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1 Introduction 



There are basically two main approaches to define the mark to market of a contingent claim: one relying on 
the no-arbitrage assumption and the other related to a hedging portfolio, those two approaches converging 
in the specific case of complete markets. A simple introduction to the different hedging and pricing models 
in incomplete markets can be found in chapter 10 of [l6j . 

The fundamental theorem of Asset Pricing [18] implies that a pricing rule without arbitrage that moreover 
satisfies some usual conditions (linearity non anticipativity . . . ) can always be written as an expectation 
under a martingale measure. In general, the resulting price is not linked with a hedging strategy except in 
some specific cases such as complete markets. More precisely, it is proved [18] that the market completeness 
is equivalent to uniqueness of the equivalent martingale measure. Hence, when the market is not complete, 
there exist several equivalent martingale measures (possibly an infinity) and one has to specify a criterion 
to select one specific pricing measure: to recover some given option prices (by calibration) [27]; to simplify 
calculus and obtain a simple process under the pricing measure; to maintain the structure of the real world 
dynamics; to minimize a distance to the objective probability (entropy [26]) ...In this framework, the 
difficulty is to understand in a practical way the impact of the choice of the martingale measure on the 
resulting prices. 

If the resulting price is in general not connected to a hedging strategy, yet it is possible to consider the 
hedging question in a second step, optimizing the hedging strategy for the given price. In this framework, 
one approach consists in deriving the hedging strategy minimizing the global quadratic hedging error under 
the pricing measure where the martingale property of the underlying highly simplifies calculations. This 
approach, is developed in [16], in the case of exponential-Levy models: the optimal quadratic hedge is 
then expressed as a solution of an integro-differential equation involving the Levy measure. Unfortunately, 
minimizing the quadratic hedging error under the pricing measure has no clear interpretation since the 
resulting hedging strategy can lead to huge quadratic error under the objective measure. 
Alternatively, one can define option prices as a by product of the hedging strategy. In the case of complete 
markets, any option can be replicated perfectly by a self-financed hedging portfolio continuously rebalanced, 
then the option hedging value can be defined as the cost of the hedging strategy. When the market is not 
complete, it is not possible, in general, to hedge perfectly an option. One has to specify a risk criteria, and 
consider the hedging strategy that minimizes the distance (in terms of the given criteria) between the pay-off 
of the option and the terminal value of the hedging portfolio. Then, the price of the option is related to the 
cost of this imperfect hedging strategy to which is added in practice another prime related to the residual 
risk induced by incompleteness. 

Several criteria can be adopted. The aim of super-hedging is to hedge all cases. This approach yields in 
general prices that are too expensive to be realistic [21] . Quantile hedging modifies this approach allowing for 
a limited probability of loss [23] . Indifference Utility pricing introduced by [29] defines the price of an option 
to sell (resp. to buy) as the minimum initial value s.t. the hedging portfolio with the option sold (resp. 
bought) is equivalent (in term of utility) to the initial portfolio. Quadratic hedging is developed in [48], [50] . 
The quadratic distance between the hedging portfolio and the pay-off is minimized. Then, contrarily to the 
case of utility maximization, losses and gains are treated in a symmetric manner, which yields a fair price 
for both the buyer and the seller of the option. 

In this paper, we follow this last approach and our developments can be used in both the no-arbitrage value 
and the hedging value framework: either to derive the hedging strategy minimizing the global quadratic 
hedging error under the objective measure, for a given pricing rule; or to derive both the price and the 



2 



hedging strategy minimizing the global quadratic hedging error under the objective measure. 
We spend now some words related to the global quadratic hedging approach which is also called mean- 
variance hedging or global risk minimization. Given a square integrable r.v. H, we say that the pair (Vq, ip) 
is optimal if (c,v) — (Vb,¥>) minimizes the functional E \H — c — J Q T vdSj . The price Vq represents the 
price of the contingent claim H and <p is the optimal strategy. 

Technically speaking, the global risk minimization problem, is based on the so-called Follmer-Schweizer 
decomposition (or FS decomposition) of a square integrable random variable (representing the contingent 
claim) with respect to an (jF t )-semimartingale S = M + A modeling the asset price. M is an (.F t )-local 
martingale and A is a bounded variation process with Aq = 0. Mathematically, the FS decomposition, 
constitutes the generalisation of the martingale representation theorem (Kunita-Watanabe representation) 
when S is a Brownian motion or a martingale. Given square integrable random variable H, the problem 
consists in expressing H as Hq + f Q T ^dS + Lt where £ is predictable and Lt is the terminal value of an 
orthogonal martingale L to M, i.e. the martingale part of S. The seminal paper is (24] where the problem 
is treated in the case that S is continuous. In the general case S is said to have the structure condition 
(SC) condition if there is a predictable process a such that At = J* a s d(M) s and J Q T a 2 d(M) s < oo a.s. In 
the sequel most of contributions were produced in the multidimensional case. Here for simplicity we will 
formulate all the results in the one-dimensional case. 

An interesting connection with the theory of backward stochastic differential equations (BSDEs) in the sense 
of [39], was proposed in [48]. [39] considered BSDEs driven by Brownian motion; in [48] the Brownian motion 
is in fact replaced by M. The first author who considered a BSDE driven by a martingale was [11]. Suppose 
that Vt = J *a s <i(M) s . The BSDE problem consists in finding a triple (V,£, L) where 



and L is an (^ r t )-local martingale orthogonal to M. 

In fact, this decomposition provides the solution to the so called local risk minimization problem, see [24] . 
In this case, Vt represents the price of the contingent claim at time t and the price Vq constitutes in fact 
the expectation under the so called variance optimal measure (VOM), as it will be explained at Section [2751 
with references to [H], [3] and [2]. 

In the framework of FS decomposition, a process which plays a significant role is the so-called mean variance 
tradeoff (MVT) process K. This notion is inspired by the theory in discrete time started by [46]; in the 
continuous time case K is defined as K t — f Q a 2 .d(M) s , t G [0, T}. [48] shows the existence of the mean- 
variance hedging problem if the MVT process is deterministic. In fact, a slight more general condition was 
the (ESC) condition and the EMVT process but we will not discuss here further details. We remark that in 
the continuous case, treated by [24] ; no need of any condition on K is required. When the MVT process is 
deterministic, [48] is able to solve the global quadratic variation problem and provides an efficient relation, 
see Theorem 15.21 with the FS decomposition. He also shows that, for the obtention of the mentioned relation, 
previous condition is not far from being optimal. The next important step was done in [38] where under the 
only condition that K is uniformly bounded, the FS decomposition of any square integrable random variable 
admits existence and uniqueness and the global minimization problem admits a solution. 
More recently has appeared an incredible amount of papers in the framework of global (resp. local) risk 
minimization, so that it is impossible to list all of them and it is beyond our scope. Two significant papers 
containing a good list of references are [51], [7] and [12]. The present paper puts emphasis on processes with 
independent increments (PII) and exponential of those processes. It provides explicit FS decompositions 
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when the process S is of that type when the contingent claims are provided by some Fourier transform (resp. 
Laplace- Fourier transform) of a finite measure. Some results of |3T] concerning exponential of Levy processes 
are generalized trying to investigate some peculiar properties behind and to consider the case of PII with 
possibly non stationary increments. The motivation came from hedging problems in the electricity market. 
Because of non-storability of electricity, the hedging instrument is in that case, a forward contract with 
value S$ = e ^ r(T ^ t) (F. f d - F Td ) where F? d is the forward price given at time t < T d for delivery of lMWh 
at time Td. Hence, the dynamic of the underlying S° is directly related to the dynamic of forward prices. 
Now, forward prices show a volatility term structure that requires the use of models with non stationary 
increments and motivates the generalization of the pricing and hedging approach developed in [31] for Levy 
processes to the case of PII with possibly non stationary increments. 

The paper is organized as follows. After this introduction and some generalities about semimartingales, we 
introduce the notion of FS decomposition and describe local and global risk minimization. Then, we examine 
at Chapter 3 (resp. 4) the explicit FS decomposition for PII processes (resp. exponential of PII processes). 
Chapter 5 is devoted to the solution to the global minimization problem and Chapter 6 to the case of a model 
intervening in the electricity market. Chapter 7 is devoted to simulations. This paper will be followed by a 
companion paper, i. e. [28] which concentrates on the discrete time case leaving more space to numerical 
implementations . 

2 Generalities on semimartingales and Follmer-Schweizer decom- 
position 

In the whole paper, T > 0, will be a fixed terminal time and we will denote by {VL,T, (J r t)te[o,T], P) a filtered 
probability space, fulfilling the usual conditions. 

2.1 Generating functions 

Let X = {X t ) t £[o,T] be a real valued stochastic process. 

Definition 2.1. The characteristic function of (the law of) Xt is the continuous mapping 

<p Xt : R -» C with tp Xt (u) = E[e iuXt ] . 

In the sequel, when there will be no ambiguity on the underlying process X, we will use the shortened notation 

<Pt for ipx t ■ 

Definition 2.2. The cumulant generating function of (the law of) X t is the mapping z i— » Log(E[e zJft ]) 
where Log(w) = log(|u;|) + iArg(w) where Arg(w) is the Argument of w, chosen in ] — tt, it]; Log is the 
principal value logarithm. In particular we have 

K Xt :D^C with e Kx ^ z} = E[e zXt ] , 

where D := {z e C | ¥,[e Re ^ Xt ] < oo, Vi e [0,T]}. 

In the sequel, when there will be no ambiguity on the underlying process X , we will use the shortened 
notation Kt for nx t ■ 

We observe that D includes the imaginary axis. 
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Remark 2.3. 1. For all z £ D, Kt(z) = Kt(z) , where z denotes the conjugate complex of z £ C. Indeed, 
for any z £ D, 

exp(K t (f)) = E[cxp(zA 4 )] = E[exp(zA t )] = E[exp(zX t )} = cxp(K t (z)) = exp(«t(z)) . 
2. For all z £ DDR , n t (z) £ R . 

2.2 Semimartingales 

An (F t )-semimartingale X = (X t ) t £[o,T] is a process of the form X = M + A, where M is an (F t )-local 
martingale and A is a bounded variation adapted process vanishing at zero. ||A||t will denote the total 
variation of A on [0, T) . Given two (Ft)~ local martingales M and N, (M, N) will denote the angle bracket of 
M and N, i.e. the unique bounded variation predictable process vanishing at zero such that AdN — (M,N) 
is an (Ft)-local martingale. If X and Y are (Ft)-semimartingales, [X, Y] denotes the square bracket of X 
and Y, i.e. the quadratic covariation of X and Y. In the sequel, if there is no confusion about the underlying 
filtration (.Ft), we will simply speak about semimartingales, local martingales, martingales. All the local 
martingales admit a cadlag version. By default, when we speak about local martingales we always refer to 
their cadlag version. 

More details about previous notions are given in chapter 1.1. of |35J. 

Remark 2.4. 1. All along this paper we will consider C-valued martingales (resp. local martingales, 
semimartingales). Given two C-valued local martingales M 1 ,M 2 then Ad 1 , Ad 2 are still local martin- 
gales. Moreover (Ad 1 , Ad 2 ) = (Ad 1 , Ad 2 ) . 

2. If Ad is a C-valued martingale then (Ad, Ad) is a real valued increasing process. 

Theorem 2.5. (-Xt)te[o,Tl * s a rea l semimartingale iff the characteristic function, t i— > ipt(u), has bounded 
variation over all finite intervals, for all u el, 

Remark 2.6. According to Theorem 1.4-18 of J35\/ . any local martingale Ad admits a unique (up to indis- 
tinguishability) decomposition, 

Ad = Ad + AI C + Ad d , 

where Ad^ = Mq = 0, Ad c is a continuous local martingale and Ad d is a purely discontinuous local martingale 
in the sense that (N,Ad d ) = for all continuous local martingales N. Ad c is called the continuous part of 
Ad and AI d the purely discontinuous part. 

Definition 2.7. An (Tt)-special semimartingale is an (Ft) -semimartingale X with the decomposition 
X = Ad + A, where Ad is a local martingale and A is a bounded variation predictable process starting at zero. 

Remark 2.8. The decomposition of a special semimartingale of the form X = Ad + A is unique, see \3S§ 
definition 4.22. 

For any special semimartingale X we define 

||X||f 2 =E[[M,M] T ]+E(|L4|||) . 

The set S 2 is the set of (Ft)-special semimartingale X for which ||-X"||| 2 is finite. 

A truncation function defined on K is a bounded function h : R — > R with compact support such that 
h(x) = x in a neighbourhood of 0. 
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An important notion, in the theory of semimartingales, is the notion of characteristics, defined in definition 
II. 2. 6 of |35| . Let X = M + A be a real-valued semimartingale. A characteristic is a triplet, (b, c, v), 
depending on a fixed truncation function, where 

1. b is a predictable process with bounded variation; 

2. c = (M c , M c ), M c being the continuous part of M according to Remark l2~6l 



3. v is a predictable random measure on E + x R, namely the compensator of the random measure /i x 
associated to the jumps of X. 

2.3 Follmer-Schweizer Structure Condition 

Let X = (Xt) te [o,T] be a real-valued special semimartingale with canonical decomposition, 

X = M + A . 

For the clarity of the reader, we formulate in dimension one, the concepts appearing in the literature, see 
e.g. [48] in the multidimensional case. 

Definition 2.9. For a given local martingale M , the space L 2 (M) consists of all predictable R-valued pro- 
cesses v = (vt)te[o,T] such that 

T \v s \ 2 d{M) s 
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(i 



< oo 



For a given predictable bounded variation process A, the space L 2 (A) consists of all predictable R-valued 
processes v = (i>t)te[o,Tl such that 

( r \v s \d\\A\\ s ^ 



E 

Finally, we set 



o 



< oo 



9 := L 2 {M)C\L 2 {A) . 
For any v £ 0, the stochastic integral process 

G t {v) := [ v s dX s , for all t G [0,T] , 
Jo 

is therefore well-defined and is a semimartingale in S 2 with canonical decomposition 



G t {v) = [ v s dM s + [ v s dA s , for all t € [0,T] . 
Jo Jo 



We can view this stochastic integral process as the gain process associated with strategy v on the underlying 
process X. 



Definition 2.10. The minimization problem we aim to study is the fol 
Given H g C 2 , an admissible strategy pair (Vo,<p) will be called optimal if (c,v) = (Vb,y) minimizes the 
expected squared hedging error 

E[(H -c-G T (v)) 2 } , (2.1) 

over all admisible strategy pairs (c, v) <G M x 9. Vq will represent the price of the contingent claim H at 
time zero. 



Definition 2.11. Let X = (Xt)te\o,T\ be a real-valued special semimartingale. X is said to satisfy the 
structure condition (SC) if there is a predictable R-valued process a = (at)te[o.T] such that the following 
properties are verified. 

1. A t = f* a s d (M) s , for all t G [0, T], so that cL4 < d (M) . 

2. [ a 2 s d{M) s < oo , P-a.s. 
Jo 

Definition 2.12. From now on, we will denote by K = {Kt)te\o,T] the cadlag process 

K t = f a 2 s d(M) s , for all t G [0,T] . 
Jo 

This process will be called the mean-variance tradeoff (MVT) process. 
Remark 2.13. In 14 the process (Kt)te[o,T] * s denoted by (-^t)ie[o,Tl- 

We provide here a technical proposition which allows to make the class of integration of X explicit. 
Proposition 2.14. If X satisfies (SC) such that E[K T ] < oo, then 9 = L 2 (M). 

Proof. Assume that E[if T ] < oo, we will prove that L 2 (M) C L 2 (A). Let us consider a process v G L 2 (M), 
then 




□ 

Schweizer in [48] also introduced the extended structure condition (ESC) on X and he provided the 
Follmer-Schweizer decomposition in this more extended framework. We recall that notion (in dimension 1). 
Given a real cadlag stochastic process X, the quantity AX t will represent the jump X t — X t ~. 

Definition 2.15. Let X = (X t ) t £[o ,t] be a real-valued special semimartingale. X is said to satisfy the 
extended structure condition (ESC) if there is a predictable R-valued process a = {at)te[o.T] with the 
following properties. 

1. A t = /„' a s d (M) s , for all t G [0, T] , so that dA < (dM) . 

2. The quantity 

is finite. 

If condition (ESC) is fulfilled, then the process 

Kt'-= f j- f', n , forallte[0,T], 
is well-defined. It is called extended mean-variance tradeoff (EMVT) process. 
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Remark 2.16. 1. (SC) implies (ESC). 

2. If (M) is continuous then (ESC) and (SC) are equivalent and K = K . 

3. K t = f* J^-d (M) s = ^ IT ^dK s , for all t e [0, T] . 

4. K t = / —~rdK s , for all t e [0, T] . 

Jo 1 — Aa s 

5. If K is deterministic then K is deterministic. 

The structure condition (SC) appears quite naturally in applications to financial mathematics. In fact, 
it is mildly related to the no arbitrage condition. In fact (SC) is a natural extension of the existence of an 
equivalent martingale measure from the case where X is continuous. Next proposition will show that every 
adapted continuous process X admitting an equivalent martingale measure satisfies (SC). In our applications 
(ESC) will be equivalent to (SC) since in Section [3721 and Section [4721 (A/) will always be continuous. 

Proposition 2.17. Let X be a (P,J-t) continuous semimartingale. Suppose the existence of a locally equiv- 
alent probability Q ~ P under which X is an (Q , Tt) -local martingale, then (SC) is verified. 

Proof. Let (-Dt)*e[o,T] be the strictly positive continuous Q-local martingale such that dP = DxdQ. By 
Theorem VIII. 1.7 of (42], M = X — (X,L) is a continuous P-local martingale, where L is the continuous 
Q-local martingale associated to the density process i.e. 

D t = cxp{L t - -{L) t } , for all t G [0, T] . 
According to Lemma IV. 4. 2 in [42], there is a progressively measurable process R such that for all t £ [0, T], 

L t = f R s dX s +O t and f R 2 s d(X) s <oo, Q - a.s. , 
Jo Jo 

where O is a Q-local martingale such that (X, O) = 0. Hence, 

(X,L) t = f R s d(X) s and X t = M t + f R s d[X] s , for all t <G [0, T\. 
Jo Jo 

We end the proof by setting at = -. ' — = Rt ■ □ 

d\X)t 

2.4 Follmer-Schweizer Decomposition and variance optimal hedging 

Throughout this section, as in Section [2~3| X is supposed to be an (J-t)-special semimartingale fulfilling the 
(SC) condition. 

We recall here the definition stated in Chapter IV. 3 p. 179 of [40] . 

Definition 2.18. Two {Tt) -martingales M,N are said to be strongly orthogonal if MN is a uniformly 
integrable martingale. 

Remark 2.19. IfM,N are strongly orthogonal, then they are (weakly) orthogonal in the sence that M[MtNt] = 
. 

Lemma 2.20. Let M, N be two square integrable martingales. Then M and N are strongly orthogonal if 
and only if (M, N) = 0. 
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Proof. Let S(M) be the stable subspace generated by M. S(M) includes the space of martingales of the form 

M( := [ f(s)dM s , for all t G [0, T] , 
Jo 

where / G L 2 (dM) is deterministic. According to Lemma IV. 3. 2 of [40], it is enough to show that, for any 
/ G L 2 (dM), g G L 2 (dN), M f and N 9 are weakly orthogonal in the sense that E[M^N%,] = 0. This is clear 
since previous expectation equals 

E[(Mf,N") T ]=E[ [ fgd{M,N)\ =0 



if (M, N) = 0. This shows the converse implication. 

The direct implication follows from the fact that MN is a martingale, the definition of the angle bracket 
and uniqueness of special semimartingale decomposition. □ 

Definition 2.21. We say that a random variable H G £ 2 (fi,JF, P) admits a Follmer-Schweizer (FS) 
decomposition, if it can be written as 



H = H a + / £f dX s + L% , P - a.s. 
Jo 



(2-2) 



where Hq G K is a constant, ^ H G and L H = [Lf )te[o.T] is a square integrable martingale, with E[Lq] = 
and strongly orthogonal to M . 

We formulate for this section one basic assumption. 

Assumption 1. We assume that X satisfies (SC) and that the MVT process K is uniformly bounded in t 
and u>. 

The first result below gives the existence and the uniqueness of the Follmer-Schweizer decomposition for 
a random variable H G L 2 {Q.,T,P). The second affirms that subspaces Gt(&) and {C 2 {Tq) + Gt(B)} are 
closed subspaces of C 2 . The last one provides existence and uniqueness of the solution of the minimization 
problem (Pi We recall Theorem 3.4 of [38] , 



Theorem 2.22. Under Assumption^ every random variable H G C 2 (fl,J r ,V) admits a FS decomposition. 

Moreover, this decomposition is unique in the following sense: 

If 

H = H a + £?dX a +L$ = H Q + / £ s H dX s + L#, 
Jo Jo 

where (Hq, £ h , L H ) and (H , £ H , L H ) satisfy the conditions of the FS decomposition, then 



H 


- H , 


P-a.s. 




= t H 


in L 2 (M) 


Lrp 




P-a.s. 



We recall Theorem 4.1 of [38] . 

Theorem 2.23. Under Assumption^ the subspaces Gt{&) and {£ 2 (J-o) + Gt(@)} are closed subspaces of 
C 2 . 
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So we can project any random variable H £ C 2 on Gt(9). By Theorem 12. 22\ we have the uniqueness of 
the solution of the minimization problem l|2.ip . This is given by Theorem 4.6 of [38], which is stated below. 

Theorem 2.24. We suppose Assumption]^ 

1. For every H £ C 2 (fl,J-,P) and every c £ C 2 (J-q), there exists a unique strategy (p^ £ 9 such that 

E[(H -c- G T (<p {c) )) 2 } = minE[(H - c - G T (v)) 2 } . (2.3) 

2. For every H £ £ 2 (f2,.F, P) there exists a unique (c^,<p^) £ C 2 (Tq) x 6 such that 

E[(H C W - G T {^)) 2 ] = mm E[(H c - G T {v)) 2 ] . 

(c,ii)e£ 2 (j r ) xe 

Next theorem gives the explicit form of the optimal strategy (p( c \ which is valid even in the case where 
X satisfies the extended structure condition (ESC). For the purpose of the present work, this will not be 
useful, see considerations following Remark l2.16l 2. 

From Follmer-Schweizer decomposition follows the solution to the global minimization problem l|2.1|) . 

Theorem 2.25. Suppose that X satisifies (SC) and that the MVT process K of X is deterministic. If H £ C 2 
admits a FS decomposition of type $2.2)) . then the minimization problem \2. 3\) has a solution ip^ £ 9 for 
any c£R, such that 

= ^ + ~ x + ^ Kt (H t - - c - G t -(<pV)) , forallt£[0,T] (2.4) 
where the process (i?t)te[o.T] is defined by 

H t :=H Q + f tfdX s + Lf , (2.5) 
and the process a is the process appearing in Definition \2.11\ of (SC). 

Proof. Theorem 3 of [48] states the result under the (ESC) condition. We recall that (SC) implies (ESC), 
see Remark 12 .161 and the result follows. □ 

To obtain the solution to the minimization problem (|2.1j) . we use Corollary 10 of [48] that we recall. 

Corollary 2.26. Under the assumption of Theorem \2.25[ the solution of the minimization problem A2.1\) is 
given by the pair (Hq, Lp( H °)) . 

In the sequel of this paper we will only refer to the structure condition (SC) and to the MVT process K. 
The definition below can be found in section II. 8 p. 85 of |40J. 

Definition 2.27. The Doleans-Dade exponential of a semimartingale X is defined to be the unique 
cadlag adapted solution Y to the stochastic differential equation, 

dY t = Yt-dX t , for all t £ [0, T] with Y = 1 . 

This process is denoted by £(X). 

This solution is a semimartingale given by 

£(X) t = cxp (X t -X - \X] t /2) + AX S ) exp(-AX s ) . 

s<t 

Theorem below is stated in [48]. 
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Theorem 2.28. Under the assumptions of Theorem \2.25l for any eel, we have 

minE[(tf - c - G T {v)) 2 ] = £{-K T ) ( (H cf + E[(L^) 2 ] + f —±~-d (E[(L H )J) ] . (2.6) 
" £S \ Jo c(—K s ) ) 

Proof. See the proof of Corollary 9 of [48] with Remark l2.16l □ 
Corollary 2.29. If(M,M) is continuous 



minE[(i? - c - G T (v)) 2 } = cxp(-K T ) ((H - cf + E[(L^) 2 ]) 



+E 



/ exp{-(K T -K s )}d(L H )^ 
Jo 



(2.7) 



Proof. Remark 12.161 implies that K = K. Since K is continuous and with bounded variation, its Doleans- 



Dade exponential coincides with the classical exponential. The result follows from Theorem 12.281 

□ 

In the sequel, we will find an explicit expression of the FS decomposition for a large class of square 
integrable random variables, when the underlying process is a process with independent increments, or is 
an exponential of process with independent increments. For this, the first step will consist in verifying (SC) 
and the boundedness condition on the MVT process, see Assumption [TJ 

2.5 Link with the equivalent signed martingale measure 
2.5.1 The Variance optimal martingale (VOM) measure 

Definition 2.30. 1. A signed measure, Q, on (fl, Tt), is called a signed 6-martingale measure, if 

(a) Q(n) = 1 ; 

(b) Q<P with ^ eC 2 (P) ; 

(c) E[^G T (v)] = for all «e9. 

We denote by P s (9), the set of all such signed Q-martingale measures. Moreover, we define 

P e (0) := {Q <G P s (0) | Q ~ P and Q is a probability measure} , 
and introduce the closed convex set, 

V d :={De£ 2 (P) | D = ^- for some Q e P s (6)} . 

dP 

2. A signed martingale measure P £ P s (0) is called variance-optimal martingale (VOM) measure if 

dP 

D = argminu e x> d Var[D 2 } = argmin£, e -p d (K[D 2 ] — l), where D = . 

The space Gt(O) := {Gt(v) \ v e 0} is a linear subspace of C 2 (P). Then, we denote by Gt{Q) ± its 
orthogonal complement, that is, 

G T (6) ± := {D G C 2 {P) | E[DG T (v)} = for any v e 0} . 

Furthermore, Gt(©)" L_l denotes the orthogonal complement of Gt(0)* l , which is the £ 2 (P)-closure of Gt(0). 
A simple example when P e (0) is non empty is given by the following proposition, that anticipates some 
material treated in the next section. 
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Proposition 2.31. Let X be a process with independent increments such that 

• X t has the same law as —X tl for any t <G [0, T]; 

• | belongs to the domain D of the cumulative generating function (t,z) i— > n t (z). 
Then, there is a probability Q ~ P such that S t = exp(A t ) is a martingale. 

Proof. For all t G [0, T], we set D t = exp{— — «t(— g)}. Notice that £> is a martingale so that the measure 
Q on (CI, Ft) defined by dQ = DTdP is an (equivalent) probability to P. On the other hand, the symmetry 
of the law of X t implies for all t G [0,T], 

S t D t = exp{^ - = cx P{-y _ «t(^)} ■ 

So S-D is also a martingale. According to [35], chapter III, Proposition 3.8 a), S is a Q-martingale and so S 
is a Q-martingale. □ 

Example 2.32. Let Y be a process with independent increments. We consider two copies Y 1 ofY and Y 2 
of -Y. We set X = Y 1 + Y 2 . Then X has the same law of -X. 

For simplicity, we suppose from now that Assumption [TJ is verified, even if one could consider a more 
general framework, see |3j Therorem 1.28. This ensures that the linear space Gt(@) is closed in C 2 (Cl), 
therefore G T (@) = G T (6) = G T (&) ±A -. Moreover, Proposition [2TT41 ensures that 6 = L 2 (M). We recall an 
almost known fact cited in [3]. For completeness, we give a proof. 

Proposition 2.33. P s (9) ^ is equivalent to 1 £ G T (Q) ■ 

Proof. Let us prove the two implications. 

• Let Q G P g (G). If 1 G Gt(6), then Q(Cl) = E Q (1) = which leads to a contradiction since Q is a 
probability. Hence 1 ^ Gt(6). 

• Suppose that 1 ^ Gt(O). We denote by / the orthogonal projection of 1 on Gr(@). Since E[/(l — /)] = 
0, then E[l - /] = E[(l - ,/) 2 ]. Recall that 1 ^ / G G T (6), hence we have E[/] 1. Therefore, we can 
define the signed measure P by setting 

P(A) = J^DdP , with D = T^-4j] ■ ( 2 - 8 ) 

We check now that P G P s (0). 

- Trivially P(Cl) = 1(5) = 1 ; 

- P <C P , by construction. 

- Let v G 6, E[DG T (v)} = t _ X E (E[(l - f)G T (v)]) = , since 1 - / G Gt^. 
Hence, P G P s (0) which concludes the proof of the Proposition. 



□ 



Remark 2.34. If 1 is orthogonal to G T (6), then f = and P G P s (6) so F S (Q) ^ 0. 

In fact, P constructed in the proof of Proposition 12.331 coincides with the VOM measure. 
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Proposition 2.35. Let P be the signed measure defined in (2.8\) . Then, 

D = argminE[.D 2 ] = argmin Var[.D] . 

D£T> d D£V d 

Proof. Let D e V d and Q such that dQ = DdP. We have to show that E[D 2 } > E[D 2 }. We write 
E[D 2 } = E[(D - Df] + E[D 2 } + ^-A— E[(J> - D)(l - /)] . 

Moreover, since / € Gt(6) yields 

E[(D - 5)(1 - /)] = E[D] - E[D] - E[Df] + E[Df] , 
= Q(Q)-Q(0). 
= . 



□ 



Remark 2.36. 1. Arai JBj gives sufficient conditions under which the VOM measure is a probability, see 
Theorem 3.4 in JE/. 

2. Taking in account Proposition 1 2. 331 the property 1 ^ Gt(&) may be viewed as non-arbitrage condition. 
In fact, in J18\/ . the existence of a martingale measure which is a probability is equivalent to a no free 
lunch condition. 

Next proposition can be easily deduced for a more general formulation, see [51] . 

Proposition 2.37. We assume Assumption^ Let H S C 2 (Vl) and consider the solution (c H ,(p H ) of the 
minimization problem i2.1\) . Then, the price c H equals the expectation under the VOM measure P of H. 

Proof. We have 

H = c H +G T ((p H )+R , 

where R is orthogonal to Gt(©) and E[R] = 0. Since P <G P s (0), taking the expectation with respect to P, 
denoted by E we obtain 

E[H] = c H + E[R\ . 

From the proof of Proposition I2.33[ we have 

~ = E t( 1 Z = 1 r E [R] - E[fR}) . 

1 1 1 -E[/] 1 -E[/] v L J u u 

Since / G G T (6) and R is orthogonal to G T (©), we get E[R] = . □ 



3 Processes with independent increments (PII) 

This section deals with the case of Processes with Independent Increments. The preliminary part recalls 
some useful properties of such processes. Then, we obtain a sufficient condition on the characteristic function 
for the existence of the FS decomposition. Moreover, an explicit FS decomposition is derived. 
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Beyond its own theoretical interest, this work is motivated by its possible application to hedging and 
pricing energy derivatives and specifically electricity derivatives. Indeed, one way of modeling electricity 
forward prices is to use arithmetic models such as the Bachelier model which was developed for standard 
financial assets. The reason for using arithmetic models, is that the usual hedging intrument available on 
electricity markets are swap contracts which give a fixed price for the delivery of electricity over a contracted 
time period. Hence, electricty swaps can be viewed as a strip of forwards for each hour of the delivery period. 
In this framework, arithmetic models have the significant advantage to yield closed pricing formula for swaps 
which is not the case of geometric models. 

However, in whole generality, an arithmetic model allows negative prices which could be underisable. 
Nevertheless, in the electricity market, negative prices may occur because it can be more expensive for a 
producer to switch off some generators than to pay someone to consume the resulting excess of production. 
Still, in [6], is introduced a class of arithmetic models where the positivity of spot prices is ensured, using a 
specific choice of increasing Levy process. The parameters estimation of this kind of model is studied in [37] . 

3.1 Preliminaries 

Definition 3.1. X = (Xt)te[o,T] is a (real) process with independent increments (PII) iff 

1. X is adapted to the filtration (^t)te[o.T] and has cadlag paths. 

2. X = 0. 

3. X t — X s is independent of T s for < s <t <T. 

Moreover we will also suppose 

4- X is continuous in probability, i.e. X has no fixed time of discontinuties. 

We recall Theorem H.4.15 of [35]. 

Theorem 3.2. Let (X t ) t ^[o,T] be a real-valued special semimartingale, with Xq = 0. Then, X is a process 
with independent increments, iff there is a version (b, c, v) of its characteristics that is deterministic. 

Remark 3.3. In particular, v is a (deterministic non-negative) measure on the Borel a-field of [0,T] x M. 

From now on, given two reals a, 6, we denote by a V & (resp. a A b) the maximum (resp. minimum) 
between a and b. 

Proposition 3.4. Suppose X is a semimartingale with independent increments with characteristics (b,c,v), 
then there exists an increasing function 1 1— > a t such that 



db t <§C da 



•t > 



dc t <C da t and v{dt 1 dx) = F t (dx)da t 



(3.1) 



where F t (dx) is a non-negative kernel from ([0, T], B([0, T])) into (K,#) verifying 




(3.2) 



and 




(3.3) 
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Proof. The existence of (at) as a process fulfilling ()3.3p and F fulfilling (|3.2|) is provided by the statement 
and the proof of Proposition II. 2.9 of [35]. (|3.3|) and Theorem 13.21 guarantee that (at) is deterministic. □ 

Remark 3.5. In particular, (b t ), (c t ) and t t— > Jj Q o xB (|^| 2 A l)f(ds, dai) /*as bounded variation for any 
BeB. 

The proposition below provides the so called Levy-Khinchine Decomposition. 

Proposition 3.6. Assume that (Xt)te[o,T] is a process with independent increments. Then 

tp t ( u ) = e *t(») , for all u€R, (3.4) 
^t, is given by the Levy-Khinchine decomposition of the process X, 

u 2 r 

* t («) = iub t - —c t + J (e mx - 1 - iuh(x))F t (dx) , for allueR , (3.5) 

where B i— > F t (B) is the positive measure v([0,t] x B) which integrates 1 A |cc| 2 /or an?/ 1 G [0,T]. 

We introduce here a simplifying hypothesis for this section. 
Assumption 2. For any t > 0, X t is never deterministic. 
Remark 3.7. We suppose Assumption^ 

1. Up to a 2ni addition of n t (e), we can write ^t{u) = Kt(iu), Vu G K. From now on we will always 
make use of this modification. 

2. ft{u) is never a negative number. Otherwise, there would be u G M*, t > such that E(cos(uXt)) = — 1. 
Since cos(uXt) + 1 > a.s. then cos(uXt) = —1 a.s. and this is not possible since Xt is non- 
deterministic. 

3. Previous point implies that all the differentiability properties of u i— > <pt(u) are equivalent to those of 
u i — ► ^i(w). 

4. IfE,[\X t \ 2 } < 00, then for all iiel, *t(w) and ^ (w) exist. 

We come back to the cumulant generating function k and its domain D. 
Remark 3.8. In the case where the underlying process is a PII, then 

D:={zeC\ E[e Re(z ) x <] < 00, Vf G [0,T]} = {z G C | E[e fle(z)Xr ] < 00} . 
In fact, for given t G [0, T], 7 G K we have 

E{e~< XT ) =E(e~< Xt )E(e~< ( - XT - Xt '>) < 00. 
Since each factor is positive, and if the left-hand side is finite, then E(e' yXt ) is also finite. 

We need now a result which extends the Levy-Khinchine decomposition to the cumulant generating 
function. Similarly to Theorem 25.17 of [45] we have. 

Proposition 3.9. Let Dq = jc G M | Jj t]x{\x\>i} eCX]y (dt, dx) < cel. Then, 
1. D is convex and contains the origin. 
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2. D = DDR. 

3. If z £ C such that Re(z) £ Dq, i.e. z £ D, then 



K t (z) = zb t + —c t + / (e zx - 1 - zh{x))v(ds, dx) . (3.6) 



[0,t]xl 



Proof. 1. is a consequence of Holder inequality similarly as i) in Theorem 25.17 of [45] . 

2. The characteristic function of the law of X t is given by 1|3.5|) . According to Theorem II. 8.1 (iii) 
of Sato [45], there is an infinitely divisible distribution with characteristics (b t ,ct,F t (dx)), fulfilling 
-f*({0}) = an d /(I A x 2 )F t {dx) < oo and q > 0. By uniqueness of the characteristic function, that 
law is precisely the law of X t . By Corollary II. 11. 6, in [45], there is a Levy process (X*,0 < s < 1) 
such that L* and X t are identically distributed. We define 



i)l < <*>} 



C f = {c e K I / e cx F t {dx) < oo} and C* = {2 e C | E [cxp(i?e(zL 
J{b|>i} 

RemarklOlsays that C T = D, moreover clearly = D . Theorem V.25.17 of @5] implies D Q = DDR, 
i.e. point 2. is established. 

3. Let t £ [0,T] be fixed; let w £ D. We apply point (iii) of Theorem V.25.17 of [45] to the Levy process 
L*. 

□ 

Proposition 3.10. Let X be a semimartingale with independent increments. For all z £ D, 1 1— ► Kt{z) has 
bounded variation and 

K d t(z) < da t . (3.7) 

Proof. Using l|3.6p . it remains to prove that 



1 1 — ► / (e za: — 1 — zh(x))v(ds, dx) 

J[0,T) xl 

is absolutely continuous with respect to (da t ). We can conclude 

K t{z)= / -r L da s + — / — -da s + / da s / (e zx - 1 - zh(x)) F s (dx) , 
J da s 2 J da s J Q J R 

if we show that 

da s { \e zx - 1 - zh(x)\F s (dx) < 00 . (3.8) 



Without restriction of generality we can suppose h(x) = xh x \<i. I|3.8p can be bounded by the sum h + h + h 
where 

h= f da s [ \e zx \F s (dx) 7 I 2 = [ da s [ F s (dx) , and I 3 =[ da s [ \e zx -l-zx\F s {dx) . 

JO J\x\>l JO J\x\>l JO ^|x|<l 

Using Proposition 13.41 we have 

h = f da s [ \e zx \F s (dx) 

JO J\x\>l 



( T da. ( \e Re(z)x \F s (dx) 

Jo J\x\>l 



[0,T]x|a;|>l 



|x|>l 

\e Re( - z)x \p(ds,dx); 
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this quantity is finite because Re(z) £ Do taking into account Proposition 13.91 Concerning h we have 

h = f da s [ F s (dx) 
Jo J\x\>i 

= [ da s [ {l/\\x 2 \)F s (dx) 

JO J\x\>l 
< CLT 

because of (|3.2j) . As far as ^3 is concerned, we have 

h < e Re ^ Z - [ da s (x 2 Al)F s (dx) 

Z J[0,T]x\x\<l 

2 

again because of (|3.2|) , This concludes the proof the Proposition. 

□ 

The converse of the first part of previous corollary also holds. For this purpose we formulate first a simple 
remark. 

Remark 3.11. For every z £ D, (cxp(zX t — K t (z))) is a martingale. In fact, for all < s < t < T . we 
have 

E[cxp(z(X t - X s ))\ = exp(K t (z) - k s (z)) . (3.9) 

Proposition 3.12. Let X be a PII. Let z £ D(~)M*. (X t ) t e[o,T] * s a semimartingale iff t 1— > Kt(z) has 
bounded variation. 

Proof. It remains to prove the converse implication. 

If t 1 — > nt(z) has bounded variation then 1 1— * e Kt ^) has the same property. Remark T3. Ill savs that e zXt = 
M t e Kt ^ where (M t ) is a martingale. Finally, (e zXt ) is a semimartingale and taking the logarithm (zX t ) has 
the same property. □ 

Remark 3.13. Let z £ D. If (Xt) is a semimartingale with independent increments then (e zXt ) is neces- 
sarily a special semimartingale since it is the product of a martingale and a bounded variation continuous 
deterministic function, by use of integration by parts. 

Lemma 3.14. Suppose that {X t ) is a semimartingale with independent increments. Then for every z £ 
Int(D), 1 1— > K t (z) is continuous. 

Remark 3.15. The conclusion remains true for any process which is continuous in probability, whenever 
1 1 ► n t {z) is (locally) bounded. 

Proof of Lemma \3.14\ Since z £ Int(D), there is 7 > 1 such that jz £ D; so 

E[cxp(zjX t )] = cxp(Kt(7z)) < exp(sup(K t (7z))) , 

t<T 

because 1 1— > K t (jz) is bounded, being of bounded variation. This implies that (exp(zX t )) tS [ T ] is uniformly 
integrable. Since (X t ) is continuous in probability, then (exp(zX t )) is continuous in C . The result easily 
follows. □ 
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Proposition 3.16. The function (t,z) i— > Kt(z) is continuous. In particular, (t,z) i— ► Kt(z), t G [0,T], z 
belonging to a compact real subset, is bounded. 

Proof. • Proposition 13,91 implies that z i— > Kt(^) is continuous uniformly with respect to t G [0, T]. 
• By Lemma \3.14\ for z £ IntD, i i— ► Kt(z) is continuous. 



• To conclude it is enough to show that 1 1— ► K t (z) is continuous for every z & D. Since 13 = IntD, there 
is a sequence (z n ) in the interior of D converging to z. Since a uniform limit of continuous functions 
on [0, T] converges to a continuous function, the result follows. 

□ 

3.2 Structure condition for PII (which are semimartingales) 

Let X = (Xt)te[o,T] be a real-valued semimartingale with independent increments and Xq = 0. We assume 
that E[|X t | 2 ] < oo. We denote by <pt(u) = E[cxp(mJf t )] the characteristic function of X t and by u ^ *£t(u) 
its log-characteristic function introduced in Proposition 13.61 We recall that ift(u) = cxp(\&t(u)). 
X has the property of independent increments; therefore 

exp(raX t )/E[exp(iwX 4 )] = exp(mX t )/cxp(* t (u)) , (3.10) 

is a martingale. 

Remark 3.17. Notice that the two first order moments of X are related to the log-characterisctic function 
of X , as follows 

E[X t ] = -i*i(0) , E[X t - X s ] = -i(*i(0) - <(0)), (3.11) 
Var(X t ) = -<(0), Far(X t -X s ) = -[<(0)-<(0)] . (3.12) 

Proposition 3.18. Let X = (X*)te[o,T] be a real-valued semimartingale with independent increments. 

1. X is a special semimartingale with decomposition X = M + A with the following properties: 

(M) t = ~%(0) and A t = -i%(0) . (3.13) 
In particular 1 1— ► — ^ t (0) is increasing and therefore of bounded variation. 

2. X satisfies condition (SC) of Definition \2.1l\ if and only if 



*;(0)«<(0) and [ 
Jo 



£(o) 



|dtf a (0)|<oo. (3.14) 



7n i/iai case 



A t = a s d(M) s with a t = i j^fM: for all t G [0,T] . (3.15) 

3. Under condition {3.1$ , FS decomposition exists (and it is unique) for every square integrable random 
variable. 

In the sequel, we will provide an explicit decomposition for a class of contingent claims, under condi- 
tion |3TT4j) . 
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Proof. 1. Let us first determine A and M in terms of the log-characteristic function of X. Using (|3.1ip 
of Remark 13.171 we get 

E[X t \F.] = E[X t - X s + X s | F a \ , = E[X t -X s )+X s , 
= -i%(0) +i*' s (0) +X S , then, 
E[X t + i%(0)\T s ] = X s +W' s (0). 

Hence, (X t + i~$f t (0)) is a martingale and the canonical decomposition of X follows 

Xt = Xt + i*[(p)-M' t (p) , 

M t At 

where M is a local martingale and A is a locally bounded variation process thanks to the semimartingale 
property of X. Let us now determine (M), in terms of the log-characteristic function of X. 

M 2 = [X t +M' t (0)} 2 , 

E[M 2 \T S ] = E{(X t + i*' t (0)) 2 \F s ] , 

= E[(X S + i<(0) +X t - X s + M' t (0) - i*l(0)) 2 |^ ; 

= E[(M S +X t - X s + i(*;(0) - <(0))) 2 |^ s ] , 



sj - 



Using gni$ and (j3~T2j) of Remark GU2l yields 

M t 2 = E[(M S - E[X t - X s ] +X t - X s ) 2 } , 

= M 2 + Var(X t - X.) = M 2 S - < (0) + < (0) 

Hence, (M 2 + ^S> t (0)) is a (Tt )-martingale, and point 1. is established. 



A t = / a s d (M) a with a, = f ^Sj^ for all t e [0, T] 
Jo «t*t ( u ) 



2. is a consequence of point 1. and of Definition 12.111 

rT ( dt* 



3. follows from Theorem l2.22l In fact K t = — I — — fr(Q) . (0) is deterministic and so Assumption 

Jo \dt* a w y 

fflis fulfilled. 

□ 

3.3 Examples 

3.3.1 A continuous process example 

Let ip : [0, T] — > R be a continuous strictly increasing function, 7 : [0, T] — * R be a bounded variation function 
such that dy <C dip. We set = W^i t ) + 7(f), where W is the standard Brownian motion on R. Clearly, 
X t = Mt + 7(<), where M t — W^u)> defines a continuous martingale, such that (M) t = [M] t = V^(t). Since 
X t ~ A^(7(t), ^(t)) for all it € R and t £ [0, T], we have 

T / n , \ u 2 ib(t) 
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which yields 

%(0)=i~/(t) and %(0) = -ip(t), 

Therefore, if — — € C 2 (dip), then X satisfies condition (SC) of Definition 12.111 with 

dip 



A t = [ a s d(M) s and a t = p- 
Jo dip 



for all t € [0,T] 



3.3.2 Processes with independent and stationary increments (Levy processes) 

Definition 3.19. X = {X t )t£[o,T} is called Levy process or process with stationary and independent 
increments if the following properties hold. 

1. X is adapted to the filtration (Tt)te[o.T] and has cadlag trajectories. 

2. X = 0. 

3. The distribution of X t — X s depends only on t — s for < s < t < T. 
4- Xt — X s is independent of T s for < s < t < T . 

5. X is continuous in probability. 

For details on Levy processes, we refer the reader to [40], [45] and [35] . 
Let X = (X t )t£[o,T] be a real-valued Levy process, with X — 0. We assume that E[|X t | 2 ] < oo and we do 
not consider the trivial case where L\ is deterministic. 

Remark 3.20. 1. Since X = {Xt)te\o,T] is a Levy process then ^t(u) = t^\{u). In the sequel, we will 
use the shortened notation ^> := 

2. ^ is a function of class C 2 and \& (0) = Var(X±) which is strictly positive if X has no stationary 
increments. 

3.4 Cumulative and characteristic functionals in some particular cases 

We recall some cumulant and log-caracteristic functions of some typical Levy processes. 

Remark 3.21. 1. Poisson Case: If X is a Poisson process with intensity X, we have that n A (z) = X(e z — 
1). Moreover, in this case the set D = C. 
Concerning the log-characteristic function we have 

= A(e m - 1) , *'(0)=«A and (0) = -A, u G M. 

2. NIG Case: This process was introduced by Barndorff- Niels en in J4$- Then X is a Levy process with 
Xi ~ NIG(a, 0,5, n), with a > \(3\ > 0, S > and \i e M. We have k a (z) = \iz + S(jq — j z ) and 
72 = y/a 2 - (P + z) 2 , D =]—a — j3,a — f3[+iR . 
Therefore 



= fiiu + 8(^0- 7i„) , where ~j iu = y/ a 2 ~ (f3 + iu) 2 
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By derivation, one gets 



*'(0) = ifi + S— and *"(0) = -S(— + 

7o 7o 7o 



Which melds a = - jCWj±jg) 

W men yields a i^,^ + p) ' 

3. Variance Gamma case: Let a, f3 > 0, 8 ^ 0. //X is a Variance Gamma process with X\ ~ X^G(a, /3, 5, //) 



w«ift k (z) = uz + 5Log - — T , where Log is again the principal value complex logarithm de- 

ya-/3z-^- ) 

fined in Section 2. The expression of k a (z) can be found in \31\ or also JJEj, table IV. 4-5 in 
the particular case fx = 0. In particular an easy calculation shows that we need z e C such that 
Re(z) e] — [3 — y/ (3 2 + 2a, —0 + \J P 2 + 2a[ so that k a (z) is well defined so that 



D=]-/3- VP 2 + 2a, -p + V 'P 2 + 2a[+il 

Finally we obtain 

^(u) = fjdu + 5 Log 

After derivation it follows 



a — (3iu + \ 



*'(O)=i(0-<*/3), *"(0) = -(a 2 -/? 2 ). 

a 

3.5 Structure condition in the Levy case 

By application of Proposition 13. 181 and Remark [3.201 we get the following result. 

Corollary 3.22. Let X = M + A be the canonical decomposition of X , then for all t £ [0, T], 

(M) t = -m"(0) and A t = -ittf'(O) . (3.16) 
Moreover X satisfies condition (SC) of Definition \2.ll\ with 

r* $>'(o) 

A t = j ad{M) s with a = ^-p^y for all t 6 [0, T] . (3.17) 

Hence, FS decomposition exists for every square integrable random variable. 
Remark 3.23. We have the following in previous three examples of subsubsection\3 '.. 
1. Poisson case: a = 1. 



2. NIG process: a 

3. VG process: a = 



5(ll + P) 
[j, — 6(3 a 



a 2 -P 2 5 
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3.5.1 Wiener integrals of Levy processes 

We take X t = J*j s dA s , where A is a square integrable Levy process as in Section [3.3.21 Then, j s dA s is 
well-defined for at least 7 € £°°([0, T}). It is then possible to calculate the characteristic function and the 
cumulative function of J Q j s dA s . Let (i, z) 1— ► t^l a(z), (resp. (t, z) 1— > tn A {z)) denoting the log-characteristic 
function (resp. the cumulant generating function) of A. 

Lemma 3.24. Let 7 : [0,T] — > R &e a Borel bounded function. 

1. The log- characteristic function of X t is such that for all 

*x t («) = / $a(u7 s )(Is , where E[exp(iuX t )] = exp (*x t («)) ; 



Lei Z?a 6e i/ie domain related to n A in the sense of Definition \2.2l The cumulant generating 
of X t is such that for all z £ {z | Rezj t € -Da /or all t £ [0, T]}, 



n A (zj s )ds. 



Proof. We only prove 1. since 2. follows similarly. Suppose first 7 to be continuous, then \ s dA s is the 
limit in probability of $3jf=o 7*j C^-fy+i ~ where = to < ^i < ■•■ < t p — T is a subdivision of [0, T] whose 
mesh converges to zero. Using the independence of the increments, we have 



^■P{i^2ltj(A tj+1 -A t ,)} 



p-i 



]E[exp{i 7 t,(A ti+1 -A tj )}] 



3=0 



p-i 



]exp{* A (7t i )(*i+i ~tj)} , 



3=0 



= exp{^(t i+1 -^)*A(7t,)} 



This converges to cxp ^ J T \Ij A (7s)em , when the mesh of the subdivision goes to zero. 

Suppose now that 7 is only bounded and consider, using convolution, a sequence 7„ of continuous functions, 
such that 7„ — > 7 a.e. and sup tg r T ] |7 n (i)| < su Pte[o,T] It(*)I- We have proved that 



3 



cxp % 



7„(s)dA s 



exp 



*A(7»(s))ds 



(3.18) 



Now, ^a is continuous therefore bounded, so Lebesgue dominated convergence and continuity of stochastic 
integral imply statement 1. □ 

Remark 3.25. 1. Previous proof, which is left to the reader, also applies for statement 2. This statement 
in a slight different form is proved in JP^ 

2. We prefer to formulate a direct proof. In particular statement 1. holds with the same proof even if A 
has no moment condition and 7 is a continuous function with bounded variation. Stochastic integrals 
are then defined using integration by parts. 
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We suppose now that A is a Levy process such that Ai is not deterministic. In particular Var(Ai) 7^ 
and so * A + °- 
In this case 



*t(«)= / ^Aiu-fshsds and = / * A (u 7s )7 s 2 rfs 

Jo Jo 



So 



*;(0) = *1(0) A.ds and <(0) = <(0) /\ 2 ds. 
Jo Jo 

Condition (SC) is verified since d^(0) < d^O) with 

d*;'(o) * A (0) 7 * {7t#} Jo s s s 1*1(0)1 

3.6 Explicit Follmer-Schweizer decomposition in the PII case 
3.6.1 Preliminaries 

Let X = (X t ) t £[o.T] be a semimartingale with independent increments with log-characteristic function 
(t,u) 1— > ^((w). We assume that (X t ) t e[o.T] is square integrable and satisfies Assumption [2j 

Remark 3.26. 1. u 1— > is 0/ dass C 2 , /or any t <G [0,T] because Xt is square integrable. 

1 t 1-4 $ t (0) and < 1 — ► ^((0) ^flfe bounded variation because of Proposition \ 3.1PX Therefore, they are 
bounded. 

3. t 1 ► ^t(n) is continuous for every u <E R. In /aci, /irsi t ^ Xt is continuous in probability. Since 
M t = X t — ^(0) is a square integrable martingale and 1 1— > ^(O) is bounded, then the family (E(X 2 )) 
is bounded and so (X t ) is uniformly integrable. So t i— ► </?t(u) is continuous and the result follows by 
Assumption^ 

4. t 1— > ^ (0) is continuous. In fact, again it is enough to prove t 1— ► ^"(0) is continuous. This follows 
if we prove that (M t ) is continuous in C 2 . This is true because M is continuous in probability and for 
any N > 0, t <E [0,T], Chebyshev implies that 

and so the family (Mf ) is again uniformly integrable. 
We suppose the following. 
Assumption 3. 1. 1 1— > ^(u) is absolutely continuous with respect to d^ t (0). 
,2. -For euery h£M, we suppose that the following quantity 



K(u) r- 



f 


d*' t (u) 


'0 


d*;'(o) 



2 



d(-*.(0)) (3.19) 



is finite. 

Remark 3.27. Ifu = Q, the previous quantity (|3.19p is finite because of the (SC) condition. 
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We consider a contingent claim which is given as a Fourier transform of Xt, 

H = f{X T ) with f(x) = / e iux n(du) , for all (3.20) 

Jr 



for some finite signed measure fi. 
Assumption 4. 



/ 

JR 



K(u)d\fj,(u)\ < oo. 
Remark 3.28. We observe that the function 

(u,t) h-> exp(\& T (u) - ^t(u)) 
is uniformly bounded because the characteristic function is bounded. 

We will first evaluate an explicit Kunita-Watanabe decomposition of H w.r.t. the martingale part M of 
X. Later, we will finally obain the decomposition with respect to X. 

3.6.2 Explicit elementary Kunita-Watanabe decomposition 

By Propostion l3.18[ X admits the following semimartingale decomposition, X t = A t + Mt, where 

A t = -i%(0) and (M) t = -%(0). (3.21) 

Proposition 3.29. Let H = f(Xx) where f is of the form 113.20]) . We suppose that the PII X satisfies 
Assumptions ®[]J\and^ Then, H admits the decomposition 

( v t = V + /* Z s dM s + O t 
\V T = H, 

with the following properties. 
1. For all t 6 [0,T], 

Z t = i 



(3.22) 



, . d (*;(«) -*;(o)) 

/ e V ^777^ L cxp {^t(u) - *,(«)} M(d«) ; (3.23) 

Jr d^ t (0) 



5. O is a square integrable {J- \) -martingale such that (0,M) = ; 
3. H = Vt where {Vt)t£[o,T\ is the {Tt) -martingale defined by 

V t = E[H\T t ] = f e mXt exp {V T (u) - * t («)} ft(du) . (3.24) 

JR 



Remark 3.30. In particular, 
1. V = E[H] ; 



2. E 



< oo . 



Proof. We start with the case /i = S u (dx) for some u 6 R so that f(x) = e m;r . We consider the (J 7 t)- 
martingale V t = E[/(X T )|^"t] = E[e iuXT \F t ]- 
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1. Clearly V = E[e iuXT ] . 

2. We calculate explicitely V t , which gives 

V t = E[e iuXT \T t ] = e ™^E[e ltl(Xr - Xf) ] = exp(iuX t - * t (u)) exp(# x (u)) 
= V t exp($ T (u)) , 

where V t — exp(iuX t — *&t(u)) defines an (J r t )-martingale. 

3. We evaluate (V,M). 

Lemma 3.31. (V, M) t = -if* V s (*' ds (u) - < s (0)) . 

Proof. We evaluate E[^M t |.F s ]. Since V and M are (J r t)-martingales and using the property of inde- 
pendent increments we get 

E[V t M t \T s ] = my t M s \F s ] + E[V t (M t - M S )\F S ] , 

= M S V S + V s E[exp{iu(X t - X„) - (tf t (u) - #,(«))} (M t - M s )] , 

= M S V S + ^ e -(**(«)-*=(«))E[e itt ( x *- x ^(-^t - ■ 
Previous expectation gives 

E [ e m(x t -x,)( M( _ Ms yj = E[e' m ( Xt - Xs )(X t - X s )] + E[e™( Xt - x ^(^(0) - *' s (0))] , 

= -i-^-E[e m(Xt - x ^} + i(*l(0) - *l(0))E[e ra ( x '- x ^] , 

= -i e **(«)-*»(«)(^( u ) _ ^( u )) + »(^(0) - ^(Ojje*^")-*^ . 

Consequently, 

E[F t M t |^] = M S V S - - <(«)) + iK(*t(0) - <(0)) 

= M s V s -iV s (^' t (u)-%(0)-(^ s (u)-^M)) ■ 
This implies that (VtM t + iVt(^ t (u) — ^(O))^ is an (Tt) -martingale. Then by integration by parts, 

f t (*;(u)-*;(o))= [ v s (*' ds (u)-v' ds (o))+ A<(«)-<(o))dK . 

Jo Jo 
The second integral term of the right-hand side being a martingale, it follows that 



t Jo 

and so 



(V, M) t = -* \ V s (*' ds (u) - ^,(0)) . (3.25) 



o 



□ 
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4. We continue the proof of the Proposition 13.291 For given (Z t ) we have 



5. We want to identify 



This naturally leads to 



f ZdM, m\ = [ Z s -d{M) s = - [ Z s < s (0). 
Jo 1 1 Jo Jo 



Z s < s (0) = -z / V s (*' ds (u)-*' da (0)) 



Zs = 1 — d¥M — Vs - ■ (3 - 26) 

6. Finally, we obtain the general case, for general finite signed measure fx, similarly to the proof of 



Theorem 14.231 fin the sequel) in the case of exponential of PII processes. The use of Fubini's theorem 
is essential. 

□ 

Example 3.32. We take X = M = W the classical Wiener process. We have f s («) = — ! rp so that 
^ s (u) = —us and (u) = —s. So Z s = iuV s . We recall that 

( nT-S 

V s = E[exp(mM / r)|J r ;! ] = exp(raW s ) exp —u 2 — — — 



In particular, Vo = exp(— !L ^-) and so 



exp(mWr) = i J u exp(iuW s ) exp ^— u 



1 dW s + exp( — ) 



2 J v 2 

In fact that expression is classical and it can be derived from Clark- Ocone formula. 

3.6.3 Explicit Follmer-Schweizer decomposition 

We introduce a quantity which will be useful in the sequel. For t g [0, T],u € R we set 

Remark 3.33. 1. r\ is defined unambiguously since d (^ t (u) — ^(0)^ is absolutely continuous with re- 
spect to dUf't (0) . 

2. rj is well-defined, because for any m£I, 

^ J dWM) d(v;m ds{ ) 

is bounded by Cauchy-Schwarz, taking into account Assumption^ point 2. 

We are now able to evaluate the FS decomposition of H = f{Xr) where / is given by l|4.27p . 
We introduce now a supplementary hypothesis. 



26 



Assumption 5. The quantity 



sup (Re(ry(u, t)) < oo 

«Gsupp/i,tG[0,T] 



Theorem 3.34. Under the assumptions of Proposition HOPl and Assumption^ the FS decomposition of H 
is the following 

H t =H + [ £ s dX a + L f with H T = H (3.28) 



and 



H t = / H(u) t ^(du) , 
£(u)tn(du) , 



(3.29) 



where 

. d(*' t (u)-*' t (0)) 

au)t = 1 — Mm — H{u)t -' 

H{u) t = exp{rj(u,T)- V (u,t) + ^ T (u)-^ t (u)} e luXt . 
Proof. Using Fubini's theorem, we reduce the problem to show that 

H(u)t = H(u)o + / £(u) s dX s + L(u) t with H(u)t = cx.p(iuXr) 



(3.30) 



o 



for fixed a£l where L{u) is a square integrable martingale and (L(u), M) = 0, where M is the martingale 
part of the special semimartingale X. Notice that by equation (|3.30l) . 

H{u) t = H(u) + eft Tv{u - ds) V(u) t with V{u) t = exp(iuX t + *t(u) - ^t(u)) . 

Integrating by parts, gives 

H(u) t =H(u) - / efr^ u ' ds )v(u) r ri{u,dr)+ / ^ ds) dV(u) r . 



JO JO 

We denote again by Z(u) the expression provided by ()3.26[> . We recall that 

dV(u) r = Z(u) r dM r + dO(u) r = Z(u) r (dX r - dA r ) + dO(u) r , 

where A is given by (|3.21|) and O is a square integrable martingale strongly orthogonal to M (i.e. (M, O) = 
0). 

H(u) t = H(u) +L(u) t + f e^^ d ^Z{u) r dX r - f e £ Z(u) r (-M dr (0))- f ^ u ^V(u) r r,(u,dr) , 

Jo Jo Jo 

where 

L(u) t = f /* is) dO(u) r , 
Jo 

is a martingale strongly orthogonal to M. To conclude, we need to choose r\ so that 

rt 

Z(u) r e^ r,(u > ds) (-i*' dr {0)) = / e£ ^ u ' ds) V(u) r r](u,dr)) . 



o 
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This requires 



d(* r (0)) 

So we define 77 as in p. 27ft , □ 
3.6.4 The Levy case 

Let X be a square integrable Levy process, with characteristic function exp(^E'(u)i). In particular, ^ is of 
class C* 2 (R). We have 

— ifA = -ttV4 and r)(u,t) = t — v ' - , ; ^ . 
d*t(0) * (0) /v ; * (0) v ; 

We remark that Assumptions [2] is verified. Concerning Assumption 02 point 1. is trivial; point 2. is verified 



because K(u) = 



* {u) 



(—T^> (it)) . On the other hand Assumption [5] is verified if 



_ / (u)tf (0) \ 

* e J < 00 • 

Since *'(0) = tE[Ai] and *"(0) < 0, (l3~3il is fulfilled if 

E[Jfi]/m(*'(u)) > -00 . 

Assumption |4] is verified if 

, 2 
* (ti) d\fi(u)\ < 00 . 



(3.31) 



(3.32) 



(3.33) 



Example 3.35. We start with the toy model Xt = aWt + mt, a, m G R. We /law = — ^cr + imu 

so ^ (u) = —u<t 2 + im and Im(^f (u)) — m. Condition (|3.32|) is always verified and Condition (|3.33f) is 
verified if 



u 2 dfi(u) < 00 



(3.34) 



(|3.34p is verified for instance in the example of the beginning of subsection \3.7\ since J_ oc u 2 e u du < 00 
for c > 0. 



Remark 3.36. In the examples introduced in Remark \3.21\ we can show that u 
so i|3.32p and (|3.33p are always verified for the following reasons. 

1. Poisson case 

We have = i\e iu . 

2. NIG case 

_i 

We have $f'(u) = i/i + i6(J3 + iu) {a 2 - {(3 + iu) 2 ) 2 . Now 



* (u) 



is bounded and 



\V(u)\<2 [\n\ 2 + 25 ] 



(3 2 + u 2 



(a 2 ~(3 2 + u 2 ) 2 +4u 2 p 2 



Since \a\ > \(3\, u 1— > \*£'(u)\ is bounded. 
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3. Variance Gamma case 

We have tf'(u) = ipt u ~ iP ^ Clearly is again bounded. 




In conclusion, we can apply Theorem 13.341 and we obtain 



V{u) t 
H(u) t 



exp(iuX t + (T - t)V(u)) , 
exp ((T - *)#(«) + 77(11, T) - »?(«, t)) e 



,mX t 




3.7 Representation of some contingent claims by Fourier transforms 

In general, it is not possible to find a Fourier representation, of the form l|3.20p . for a given payoff function 
which is not necessarily bounded or integrable. Hence, it can be more convenient to use the bilateral 
Laplace transform that allows an extended domain of definition including non integrable functions. We refer 
to [17], [H] and more recently [20] for such characterizations of payoff functions. This will be done in the 
next section. However, to illustrate the results of this section restricted to payoff functions represented as 
classical Fourier transforms, we give here two simple examples of such representation extracted from [20] : 

1. A variant of the digital option is the so-called asset- or-nothing digital, where the option holder receives 
one unit of the asset, instead of currency, depending on wether the underlying reaches some barrier or 
not. Hence, the payoff of the asset-or-nothing digital put with barrier is 

r Di+iu 

n / \ f -m 1 It / \ I HIT />/ \ 1 



In both cases the measure fj, is finite. 

4 Follmer Schweizer decomposition for exponential of PII processes 

In this section, we consisder the case of exponential of PII corresponding to geometric models (such as the 
Black-Scholes model) much more used in finance than arithmetic models (such as the Bachelier model). 
The aim of this section is to generalize the results of [31] to the case of PII with possibly non stationary 
increments. Here again, this generalization is motivated by applications to energy derivatives where forward 
prices show a volatility term structure that requires the use of models with non stationary increments. 

4.1 A reference variance measure 

We come back to the main optimization problem which was formulated in Section [2j We assume that the 
process S is the discounted price of the non-dividend paying stock which is supposed to be of the form, 



where sq is a strictly positive constant and X is a semimartingale process with independent increments (PII) , 
in the sense of Definition 13. 1\ but not necessarily with stationary increments. 




2 



The payoff of a self quanto put option with strike K is 



/(*) 




S t = s cx-p(X t ) , for all t £ [0, T] , 
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Remark 4.1. Let 7 G M* , 

1. E[cxp(7(X t - X a ))] > 0, since X t - X s > -00 a.s. 

2. exp(7(X f — X s )) has a strictly positive variance if (X t — X s ) is non- deterministic. 
We introduce a new function that will be useful in the sequel. 

Definition 4.2. For any t G [0,T], let pt denote the complex valued function such that for all z,y G D 

Pt(z, y) = K t (z +y)- K t (z) - Kt{y) ■ (4-1) 

For all z G D, then z G D and pt(z,z) is well defined. To shorten notations pt will also denote the real 
valued function defined on D such that, 

Pt (z) = p t (z, z) = K t (2Re(z)) - 2Re(Kt(z)) . (4.2) 

Notice that the last equality results from Remark \2.3\ . 

An important technical lemma follows below. 

Lemma 4.3. Let z £ D, with z ^ 0, then, t 1— > pt{z) is strictly increasing if and only if X has no 
deterministic increments. 

Proof. It is enough to show that X has no deterministic increment if and only if for any < s < t < T, the 
following quantity is positive, 

pt(z) - p s (z) = [ Kt (2Re{z)) - K s (2Re(z))] - 2 Jte («*(«) - k s (z)) . (4.3) 

By Remark [3H11 for all z G D, we have 

exp[«t(z) - K a (z)] = E[exp(zA*)] , where A* = X t - X s . 



Applying this property and Remark 1231 1., to the exponential of the first term on the right-hand side of 
(|4~3]) yields 



cxp [K t (2Re(z)) - K s (2Re(z))] = E[exp(2iJe(z)A*)] = E[exp((z + f)A*)] 

= E[|exp(zA 4 s )| 2 ] . 

Similarly, for the exponential of the second term on the right-hand side difference of l|4.3p , one gets 

|2 



exp [2Re(Kt(z) - k s (z))] = exp (nt(z) - k s (z)) + (nt(z) - k s (z)) = |E[exp(zA*)]| 

Hence taking the exponential of pt{z) — p s (z) yields 

r / \ f , E[|exp(*A*)| 2 ] 

exp[p t (z)-p s (z)}~l = L-^m 1, 

|E[exp(zA*)]| 



|r " (z - >l I - 1 , where Ttfz) = exp(; 
|E[P s (z)]| 2 

Var [Re(rl(z))] + Var [Jm(r*(z))] 
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(4.4) 



• If X has a deterministic increment A* = X t — X s , then r^(z) is again deterministic and l|4.4p vanishes 
and hence t — > pt{z) is not strictly increasing. 

• If X has never deterministic increments, then the nominator is never zero, otherwise i?e(r'(z)), 
7ro(r*(z)) and therefore T*(z) would be deterministic. 

□ 

From now on, we will always suppose the following assumption. 
Assumption 6. 1. (X t ) has no deterministic increments. 
2. 2eD. 

Remark 4.4. 1. In particular for 7 € D, 7 7^ 0, the function 1 1— > pt{~f) is strictly increasing. 

. ., Var(exp(Al)) 

2. If z = 1, (|X7l) equals which is a mean-variance quantity. 

L - ' (E[ex P (A' s )]) 2 

We continue with a simple observation. 
Lemma 4.5. Lei I be a compact real interval included in D. 

supsupE[«S7] < 00 . 
7 e/ t<T 

Proof. Let t £ [0, T] and 7 G /, we have 

E[S7] = «2 exp{« t (7)} < max(l, s s Q upI ) exp( sup | Ct) I ) • 

t<T, 7 eJ 

since k is continuous. □ 

We state now a result that will help us to show that Kdt{z) is absolutely continuous with respect to 
fttt(l)=K«tt(2)-2Ml). 

Lemma 4.6. We consider two positive finite non-atomic Borel measures on E c R™, and ^. We suppose 
the following: 

1. n -c ^ ; 

2. ^ /or euen/ open 6aH 0/-E. 

Tften — — := h 7^ v a.e. In particular /1 and v are equivalent, 
dv 

Proof. We consider the Borel set 

B = {xe E\h(x) = 0} . 

We want to prove that v(B) = 0. So we suppose that there exists a constant c > such that v(B) = c > 
and another constant e such that < e < c. Since v is a Radon measure, there are compact subsets K £ and 

K± of E such that 

2 

K € cKc<zB and z/(B - K e ) < e , z/(5 - A' f ) < | . 
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Setting e = |, we have 



u(K e ) > | and i/(JT 4 ) > ^ 



By Urysohn lemma, there is a continuous function ^? : E — * R such that, < ip < 1 with 

o? = 1 on K e and w = on J\ ? . 

2 

Now 

/ ^(x)^(dx) > ^(^T e ) > ^ > . 

By continuity of </? there is an open set O C E with <p(a;) > for ieO, Clearly O C if ;l C -B; since O is 
relatively compact, it is a countable union of balls, and so B contains a ball I. The fact that h = on I 
implies p(I) = and this contradicts Hypothesis 2. of the statement. Hence the result follows. 

□ 

Remark 4.7. From now on, in this section, dpt = pat will denote the measure 

dp t = Pdt{l) = d(«t(2) - 2«t(l)) . (4.5) 

According to Remark \4-4\ 1., it is a positive measure which is strictly positive on each interval. This measure 
will play a fundamental role. 



Remark 4.8. 1. If E = [0,T], then point 2. of Lemma \4-6\ becomes p,(I) ^ for every open interval 
/C[0,T]. 

2. The result holds for every normal metric locally connected space E, provided v are Radon measures. 
Proposition 4.9. Under Assumption^ 

d{K t (z)) < dp t , for all z e D . (4.6) 



Proof. We apply Lemma 14761 with d\i = dp t and dv = da t . Indeed, Corollary 13.101 implies Condition 1. of 
Lemma l4~6l and Lemma l4~3l implies Condition 2. of Lemma 14.61 Therefore, dat is equivalent to dp t . □ 



Remark 4.10. Notice that this result also holds with dp t (y) instead of dpt = dp t (l), for any y <G D such 
that Re(y) ^ 0. 

4.2 On some semimartingale decompositions and covariations 

Proposition 4.11. Let y 7 z e D such that y + z, 2Re{y) 1 Re{y) + 1, 2Re{z) and Re{z) + 1 e D. Then S z is 
a special semimartingale whose canonical decomposition 5| = M(z) t + A{z)t satisfies 

A(z) t - / S z u _ Kdu (z) , (M(y),M(z)) t = f S y a t z p d u(z,y) , M(z) = s z , (4.7) 
Jo Jo 



where dp u (z) is defined by equation \4-1§ - In particular we have the following: 

1. (M(z),M) t = /^W^l) 

2. (M(z),M(z)) t = J*S^ Pau (z) . 
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Proof. For simplicity, we will only treat the case y = 1 in Q4.7J1 . i.e. statement [TJ The general case will follow 
similarly. By Remark T3.1H N{z) t := e~ Kt ^S^ is a martingale. Integration by parts yields 
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S * = e ^-) N ( z ) t = M(z) t + A(z) t with M (z) = s§, A(z) t = \ S z u _n du {z) and 
[M(z),M] t = [S z ,S] t} 



t r-t 

1 QZ Q\ I 02 IQl / Ol JQZ 



St^t SqSq / Ss-dS s / S s _dS 
h Jo 



s ' 



= SfSl - SftSg - f S z s _dM s - f S z _dA s - f S x s _dM{z) s - f Sl_dA(z) s , 
Jo Jo Jo Jo 

= S z+1 -S z+1 - f S z s _dM s - f Sl_dM(z) s - f S^K^i)- f S^kM , 
Jo Jo Jo Jo 

= M(z + l) t - [ S z _dM s - f Ss-dM(z) 3 + f (K ds (z + 1) - K ds (z) - Kdsil^t 1 . 
Jo Jo Jo 

Note that the first three terms on the right-hand side are local martingales. Since (M(z), M) t is the pre- 
dictable part of finite variation of the special semimartingale M{z)M, equation ((TJ) follows. 

□ 

Remark 4.12. Lemma \4-5\ implies that E[| (M(y), M(z)}\] < oo and so M(z) is a square integrable mar- 
tingale for any z E D such that 2Re(z), Re(z) + 1 E D. 

4.3 On the Structure Condition 

If we apply Proposition 14.111 with y = z = 1, we obtain S = M + A where M is a martingale and 

A t = ( S u -K du {l) , (4.8) 
Jo 

and 

(M,M) t = f S 2 u _{K du {2)-2K du {l))= [ S 2 u _p du . (4.9) 
Jo Jo 

At this point, the aim is to exhibit a predictable K- valued process a such that 

1. At= fia a d(M), ■ 

2. K T = Jq ' a 2 d(M) s is bounded. 

In that case, according Theorem [5221 there will exist a unique FS decomposition for any H E C 2 and so the 



minimization problem (|2,1|) will have a unique solution, by Theorem 12.251 
Proposition 4.13. Under Assumption® we have 



A t = f a 3 d(M) s , (4.10) 
Jo 
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where a is given by 

Xu ... . d.K u (l) 



= — — with X u := — , for all u G [0,T] . (4.11) 



Moreover the MVT process is given by 



K t = I ( %ii2T) dpu . (4.12) 



o \ d P 

Corollary 4.14. Under Assumption^ the structure condition (SC) is verified if and only if 

JO \ d Pu J 

In particular, {Kt) is deterministic therefore bounded. 

Proof of Proposition \4-TS[ By Proposition 14.9) dn t (l) is absolutely continuous with respect to dpt. Set- 
ting a u as in l|4.1ip . relation l|4.12p follows from Proposition I4.11[ expressing K t = L o? u d (M) u . □ 

Lemma 4.15. The space is constituted by all predictable processes v such that 

r-T 



E[ / v%S%_dp t ] < oo . 
Jo 



Proof. According to Proposition 12.141 the fact that K is bounded and S satisfies (SC), then »60 holds if 
and only if v is predictable and E[J Q T v\d (M, M) t ] < oo. Since 



(M,M) t = I S 2 s _d Ps , 



o 

the assertion follows. □ 

4.4 Explicit Follmer-Schweizer decomposition 

We denote by V the set of z € D such that 

dn u (z) 



dp u < oo. (4.13) 



dp u 

From now on, we formulate another assumption which will be in force for the whole section. 
Assumption 7. 1 6 D. 

Remark 4.16. 1. Because of Proposition U. !A — ^— - exists for every z 6 D. 

apt 

2. Assumption^ implies that K is uniformly bounded. 

The proposition below will constitute an important step for determining the FS decomposition of the 
contingent claim H — f{Sr) for a significant class of functions /, see Section [4751 

Proposition 4.17. Let z G T> with z + 1 G V and 2Re(z) G D. 

1. S'f G c 2 (n,T T ). 
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2. We suppose Assumptions^ and and we defir 



7 ( Z ,t):= d -^l,te[0,T). (4.14) 
apt 



J T \"f(z,t)\ 2 p dt < oo and 



7]{z,t) := Kt(z)- -y(z,s)n ds (l) 
Jo 

(4.15) 

Jo d Ps 

is well- defined and r)(z,-) is absolutely continuous with respect to p ds and therefore bounded. 

3. Under the same assumptions H{z) = Sf, admits a FS decomposition H(z) = H(z)q + ^(z) t dSt + 
L(z)t where 

H(z) t := e$ '^ z ' ds ^> St , (4.16) 

f(z)t := i(z,t)e& "<■'>*'') S?! 1 , (4.17) 

L(z) t := H(z) t -H(z) - f £(z) u dS u . (4.18) 

Jo 

Proof. 1 . is a consequence of Lemma 14.51 

2. 7(2, •) is square integrable because Assumption [7] and z, z + 1 e V. Moreover rj is well-defined since 



T 



\l{z,s)\ 



dp s 



o Jo 



Pds) < h(z,s)\ 2 p ds \—^\ 2 Pds . (4.19) 



3. In order to prove that l|4.16p . (|4.17p and <|4. 18(1 constitute the FS decomposition of H(z), taking into 
account Remark 12.191 we need to show that 

(a) H(z)o is JFo-measurable, 

(b) (L(z),M)=0, 

(c) e 6, 

(d) is a square integrable martingale. 

Point (a) is obvious. Partial integration and point [Tl of Proposition 14. Ill yield 

ft r t 



H{z) t = H{z) + / e^^ ds) dS z + / S z u d{e ! - ,l{z - ds) ) 
Jo Jo 

= H(z)o+ I' e^ ,l{z ' ds) dM{z) u + f e^ ri(z ^ ds) dA(z) u + I ' S*d(e^ ^ z ' ds ^) 



JO Jo 

t p t rt 



= H(z) Q + / e^^ z ' ds UM{z) u + / e^^'^dA(z) u - / ^ z - ds) S z r]{z, du) 
Jo Jo Jo 

= H(z) + f e^^ ds UM{z) u - [ efu^>Ws*r,(z,du)+ f e% ^ z ' ds ^ S z _n du (z) 
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On the other hand 



Jo 



£,(z) u dS u = / £,(z) u dM u + / £{z) u dA u , 
Jo Jo 

= f £(z) u dM u + f t(z) u S u -K du (l) , 
Jo Jo 

= [\(z) u dM u + f\(z,u)e^^ z - ds '>S z u _K du (l) . 
Jo Jo 

Hence, 

L(z) t = H{z) t - H{z) Q - f £(z) u dS u , 

Jo 

i-t r-t 

X v(z,ds) dM ^ u _ / e C^ z ' ds) S z r 1 (z,du) + / e J"« n{z ' ds) Sl_K du {z) 



t 



£{z) u dM u - / 7 (z,u)e^^-ds) s z_ Kdu{1) 
Jo 

fui(*> d *)dM(z) u - [ £{z) u dM u 



+ / e/« *(»>*>) S£_ [ Kdu (z) - r,(z, du) - j(z, u) Kdu (l)}. 
Jo 

Then, by definition of i] in (|4.15p . rj(z, du) = n du {z) — 7(2, u)n du (l) , hence, 

L(z) t = fe^^dM(z)„ - f\(z) u dM u , (4.20) 



which implies that L(z) is a local martingale. 
From point Q] of Proposition 14.111 it follows that 

(L(z),M) t = f ef^^d(M(z),M) u - f az) u d(M,M) u , 
Jo Jo 

= f e $Z ^ ds) Slt l p du {z, 1) - f Z{z) u Sl_ Pdu , 
Jo Jo 

= f ett^M S *+} Pdu ( z ,l)- f 1 (z,u)ef^^S z u t 1 p du ■ 
Jo Jo 

Consequently, 

(L(z),M) t = [ t ef^^S z u t 1 [p du (z 7 l)~l(z,u)p du ] . 
Jo 

Then by definition of 7 in l|4.14p . p d t{z, 1) = j(z 7 t)p dt , which yields, 

(L(z),M) t = 0. (4.21) 
Consequently, point (b) follows. To continue the proof of this proposition we need the lemma below. 
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Lemma 4.18. For all z G C as in Proposition \4-17\ Apt a.e. we have 



1. l(z,t) = 7OM) ; 



2. rj(z,t) = ij(z,t) . 

Proof. Using Remark 12.31 1) we observe z,z + leP. 

1. By definition of 7 in <|4. 14|) . 7(2, t)pdt = Pdt{z, 1) ■ Then, taking the complex conjugate of the integral 
from to t and using Remark [231 1 yields, 



j(z, s)p ds = / 7(2, s)p ds , 
Jo 



= Pt(z, 1) = «t(z + 1) - Kt(z) - «t(l) , 
= = K f (z + 1) - K t (f) - K t (l) = p t (z, 1) , 

7(z,s)p ds . 



2. By definition of 77 in (|4.15p . r)(z, t) = n t (z) — J Q 7(2, u)/td«(l) , so taking the complex conjugate, 



r}{z,t) = K t (z) - ~f(z,s)K ds (l) , 







K t( z ) - / 7fo s ) K ds(!) 



□ 



We continue with the proof of point 3. of Proposition 14.171 It remains to prove that L(z) is a square- 
integrable martingale for all z e D and that Re(£(z)) and Im(£(z)) are in 0. (|4.2Q[) says that 

L{z) t = I e^^ z ' da UM s {z)- (\(z) s dM s . 
Jo Jo 

By Proposition 14.111 and Lemma 14.181 it follows 



L(z) t = L{z) t , (4.22) 



hence, 



L(z),L(z)) t = (L(z),L(z)) t , 



(4.23) 





t 



az) s e f ^^ du) d(M 7 M(z)) s . 
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By Proposition 14. 1 II we have 

Using Lemma [4. 181 and expressions (|4. 14[) . and l|4.17p of j(z, s) and £{z) s , we have 

(L(z),Wj) t = {V 2 ^^^ , 



JO 
t r t 



e Sj2Re(n(z,du)) S 2Re{z) pds{z) _ / ^ 2Re (v(z 4u)) g^) ^ 



Consequently 



(L(z),L(zj) t = jf eX 2 ^>^)S 2 s R ^[p ds (z) - h(z, S )\ 2 Pds ] . (4.24) 



Then, point 2. implies 



rp rp 

\az)s\ 2 S 2 s _ Pds = [ Z(zU(z)sS 2 s-Pds 
JO 

1 (z,s)e^^ d ^SII 1 1 (z,s)e^^ d ^SIl 1 S 2 s _p ds , (4.25) 
\ 7 (z,s)\ 2 ef? 2R «^ du »S 2 » e ^p ds . 
Taking the expectation in <|4.25|) . using again point 2., <|4. 14)1 . (|4.15[) and Lemma l4~5l we obtain 



T 



II 



3 



L(z),L(z)j ^ < oo . (4.26) 

Therefore, L is a square-integrable martingale. Similarly, <|4.25|) yields that i?e(£(z)) g 6 and Im(£ y (z)) e 6. 
This concludes the proof of Proposition l4.17l □ 

4.5 FS decomposition of special contingent claims 

Now, we will proceed to the FS decomposition of more general contingent claims. We consider now options 
of the type 

H = f(S T ) with f(s) = f s z U(dz) , (4.27) 

Jc 

where LT is a (finite) complex measure in the sense of Rudin [44], Section 6.1. An integral representation of 
some basic European calls can be found later. 
We need now the new following assumption. 

Assumption 8. Let Iq = suppLT n M. We denote I := [inf Jo A 2 inf Jo, 2 sup Jo V sup Iq + 1] . 
1. Vz e suppLT, z,z + l eV. 
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2. I C D and sup rc6JU | 1 -j 



d(K t {x)) 
dpt 



< 00. 



Remark 4.19. 1. Point 2. of Assumption^ implies sup zeJ+iH ||Krf t (i?e(z))|| T < oo . 

2. Under Assumption^ H = f(Sr) is square integrable. In particular it admits an FS decomposition. 

3. Because of in Proposition ^ -9\ the Radon- Nykodim derivative at Point 2. of Assumption^ always 
exists. 

We need now to obtain upper bounds on z for the quantity (|4.26|1 . We will first need the following lemma. 
Lemma 4.20. There are positive constants 01,02,03 such that dp s a.e. 
1. 

dRe(ri{z,s)) 
sup V - < ci. 

zfzsuppn o>Ps 



For any z G suppll 



if m2 ^ dp s (z) 
\l(z,s)\ <— 



< c 2 - c 3 



dRe(r)(z,s)) 
dp s 



sup 

sEsuppII Jo 



2Re(r](z, dt)) exp( / Re(r)(z,ds))) < oo. 



Remark 4.21. According to Proposition 1 1— > Re(r](z, t)) is absolutely continuous with respect to dpt. 

Proof (of Lemma [4~20]) . 

The proof is inspired by Lemma 3.9 of [31]. According to Point 2. of Assumption [8] we denote 



en := sup 

x£l 



d(Kt{x)) 



dpt 



For 2 € suppll, £ e [0, T], we have 

r}(z,i) = K t (z) — J y(z,s)dK s (l) and r)(z,t) = Kt{z) — J 
Then, we get Re(r)(z, t)) = Re(K t {z)) — f Q Re('j(z, s))dK s (l) . We obtain 



Re(r}(z,ds)) < Re {k t (z) - n t {z)) + 
T Re(dn s {z)) 



dp s 



-dp a 



y(z,s)dK s (l) 



j(z,s)dn s (l) 



(4.28) 



(4.29) 



Since \ L{z),L(z)j is increasing, and taking into account l|4.24p . the measure, (dp s {z) — \^(z, s)| 2 e?p s ) , is 
non-negative. It follows that 



dps(z) 
dp s 



\l(z, s)\ 2 > , dp s a.e. 



(4.30) 
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Remark 4.22. By Lemma (14.3011 . in particular the density ^ - is non-negative dp s a.e 

dp s 



Consequently, 



dRe(K s (z)) dn s (2Re(z)) 
2 ; < ; , dp s a.e. 



dp s dp s 

In order to prove 1. it is enough to verify that, for some Cq > 0, 

dRe0n(z,s)) ldRe(K s (z)) 

3 < c + dp s a.e. 

ap s 2 dp s 

In fact, l|4.3ip and Assumption [8] point 2. and (|4.28[> . imply that 



dRe(r,(z,s)) 1 

1 < c + -en =: c\. 

dp s 2 



To prove (|4.32j) it is enough to show that 



Re(v(z, T) - Tj(z, t)) < c {p T - Pt ) + -Re(K T (z) - K t {z)), Vi g [0, T}. 



Again Assumption [8] point 2. implies that 

rT 



j(z,s)dK s (l) 



< C12 / \i(z,s)\dp 8 , 



where cyi 



Using (|4.30p . and Assumption [8] it follows 



2 . rf/Q s (z) dn(2Re(z)) 2dRe(K s (z)) 



dp s 



dp s 



dp s 



< en 



2dRe{n s {z)) 
dp s 



This implies that 



where C13 > is chosen such that cf 3 > 4cf 2 + c\ 2 c\i. Consequently 



7(z,s)dK s (l) 



dRe(K s (z)) 



dp s 



Coming back to (|4.29|) . we obtain 

Re(r](z,T) - r)(z,t)) < 



Re(dn s (z)) 
dp s 



C13 + ^ 



Re(dK s {z)) 



dp s 



< 



lRe(dn s (z)) ( Re(dn s (z)) 



dp s 



\ dp s 



dp s 

Cl3 



(|4.3ip and Assumption [8] allow to establish 



Re(r)(z,T) - r)(z,t)) < J dp s ( c 



1 dRe(K s (z)) 

2 dp s 
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where Co 



ci 3 . This concludes the proof of point 1. 



In order to prove point 2. we first observe that l|4.32p implies 



dRe(n s (z)) ( dRe(rj(z,s)) 
< 2 Co - 



dp s a.e. I|4.36p implies 



dp s 



\y{z,s)\ < c 2 i 



dp s 

L rffle(77(z, S )) 
dp a 



(4.38) 



(4.39) 



where C21 = en + 4co. Point 2. is now established with C2 = C21 and C3 = 4. 
We continue with the proof of point 3. We decompose 



where 



A+(z,t) = 



Re(r](z,t)) = A + (z,i) - A~(z,t) , 

* (dRe(r)(z,s)Y 
V dp s 



dp s , and A (z,t) = / 
Jo 



( dRe(j](z, s)) 
V dp s 



dp s 



A + (z, .) and A (z, .) are increasing non negative functions. Moreover point 1. implies 

A+(z,t) < c lPt . 



At this points for z 6 I + il 

l-T 



Re(ri(z,dt))eft ^v{*,d»)) 







< 



e 



(A~{z, dt) - A+{z, dt)) e 2 ^ *«(«»(*.*•)) 
J 4-( z ,di) e 2 ( A+ ( j; > T )-^ + ( z '*))e- 2 (^"( z ' T )- A "( ;s '* ) ) 
pt [ T e - 2 ( A ~^- A ~^A-(z,dt) 



< 



2 

g2cip T 



1 -e 



-2A~(z,T) 



} 



which concludes the proof of point 3 of Lemma 14.201 



Let 7 = sup zgJ (2Re(z)) ) by Lemma [475} it follows 



sup E 

z£l,s<T 



g2Re(z) 



< 00 . 



□ 



(4.40) 



Theorem 4.23. Let H be a finite complex-valued Borel measure on C. 

Suppose Assumptions®^^ Any complex-valued contingent claim H — /(St), where f is of the form \4.27\ ), 
and H S C 2 , admits a unique FS decomposition H = Hq + J Q T £tdSt + Lt with the following properties. 

1. H eC 2 and 

• H t = J H(z) t U(dz) , 

. Ct = H(z)Mdz), 
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L t = f L(z) t n(dz) , 



where for z G supp(II), H(z),^(z) and L(z) are the same as those introduced in Proposition and 
we convene that they vanish if z ^ supp(II). 



2. Previous decomposition is real-valued if f is real-valued. 
Remark 4.24. Taking II = 5 Zo (dz), z e C, Assumption^ is 

Proof, a) /(St) € C 2 since by Jensen, 



to the assumptions of Proposition 



E 



n{dz)s^ 



< / |n|(cfe)£;|^ ez ||n|(c) < su P £(s| x )|n|(c) 2 



where |II| denotes the total variation of the finite measure II. Previous quantity is bounded because 
of Lemma 14.171 

We go on with the FS decomposition. We would like to prove first that H and L are well defined 
square-integrable processes and E(J^ \^ s \ 2 d(M) s ) < oo. 
We denote K = supp{Yl). By Jensen's inequality, we have 



E 



L(z) t IL(dz) 



<E / \U\(dz)\L t (z)\ 2 t )\U(C)\= / |n|(dz)E[|L t (z)|?]|II| 



Similar calculations allow to show that 

E[£ 2 ] < |n|(C) J |n|dz)E[|6(z)| 2 ] and E[L 2 t ] < |n(C)| 
We will show now that 



• (Al) 

• (A2) 

• (A3) 



suPt<T,zei+iR E l\ H t(z)\ 2 } < oo : 
J c \Tl\(dz)E[\L t (z)\ 2 ]<^ 



I J dp t S 2 I \t t (z)\ 2 \ll\(dz) ) <oo. 



(Al): Since H(z) t = e$ "('-^Sf* , we have 

\H(z) t \ 2 = H{z) t llWt = ^ 2Re(v{z4s)) S 2Re(z) 



so 



E[|F(z) t | 2 ] = e£ 2R <^ z > ds »E[S 2 t Re{z) ] < e J?**><nl*.*>)) sup E[S7] , 



t<T 



with 7 = sup z(S/ 2Re(z). Inequality (|4.40[) and Lemma l4.20l imply (Al). Therefore (H t ) is a well-defined 
square-integrable process. 

(A2): E[|L t (z)| 2 ] < E[|L T (z)| 2 ] = E[^L(z),L(z)j ] , where the first inequality is due to the fact that 
|£t(z)| 2 is a submartingale. 

f-T 



E 



L(z),L(z) 



E 
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By Fubini, Lemma l4~5l and <|4.24)1 . we have 



E 



L(z),L(z) 



< 



< C4 



jf 2Re(r)(z,du) 



jJ 2Re(ri(z,du) 



[St 

dps(z) _ 
dp s 

dp s (z) 



dp s 



' dp s (z) 
dp s 



dp s , 



dp s 



E[S 2 * e{z> }d Ps 



where C4 = sup s<T EfS's^ 2 ' 1 ] . 

According to Lemma 14.201 point 2, previous expression is bounded by 04/(2), where 



I(z) 



J dptcxp^j^ 2Re(ij(z,ds)) 
C2h(z) + c 3 h(z) , 



C2 - c 3 - 



dRe{r](z,i)) 



dp t 



(4.41) 



where 



h{z) — J dptcxp^j^ 2Re{rj{z, ds))^j and Is(z) = J exp iyj^ 



2Re(r](z, dsj) dRe(r](z,ds)) 



Using Lemma T4.201 we obtain 



sup |/i(z)| < pr exp (2cipt) and sup |i2(z)| < 00 

zel+iW z£l 



and so 



sup 



L(z),L(z) 



< 00 



(4.42) 



(4.43) 



This concludes (A2). 

We verify now the validity of (A3). This requires to control 



Pdt sf / \n\{dz)Uz) t \ 



< E 



j(z, t) exp / Re(r](z, ds)) S 1 



2-1 



Using Jensen inequality, this is smaller or equal than 

i-T 



|n(c)| / \n\(dz) / Pdt E 

Jc Jo 
Lemma 14.201 gives the upper bound 



5' 



2Re{z) 



|7(z,i)| 2 cxp 2 / Re(r](z,ds)) 



|n|(c) sup ] 

t<T,~/£l 



s. 



2Re(z) 



\U\(dz)I(z) , 



where I{z) was defined in l|4.42p . Since II is finite and because of l|4.43p . (A3) is now established. 

We show now that (L t ) is an (.F t )-martmgale. Let 0<s<i<T,B6 T s . By Proposition 14.171 since 
(L[z)t) is a martingale, we obtain 



E[(L t - L s )l B } = E[ / (L(z) t - L(z) s )U(dz)l B ] 
Jc 
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By Fubini's theorem we conclude that 



E[(L t - L S )1 B ] = f E[(L(z) t - L(z) s )l B ]U(dz) 



and E[(L(z) t - L(z) s )l B ] = 0. So 



E[(L t - L S )1 B ] = . 



Hence, L is a square-integrable martingale. 

Similarly, it can be shown that E[(M t L t — M S L S )1 B ) = and so ML is a square-integrable martingale 
as well. Hence L is orthogonal to M. By Fubini's theorem for stochastic integrals, cf. [40], Theorem 
IV. 46, we have 



Consequently, (-ffo,£,L) provide a (possibly complexe) FS decomposition of H. 

b) It remains to prove that the decomposition is real-valued. Let (Ho,£,L) and (Ho,£,L) be two FS 
decomposition of H. Consequently, since H and (St) are real-valued, we have 



which implies that = Im(H ) + L Im(£ s )dS s + Im^Lr)- By Theorem I2.22[ the uniqueness of 
the real- valued Follmer-Schweizer decomposition yields that the processes (H t ),((,t) and (L t ) are real- 
valued. 



4.6 Representation of some typical contingent claims 

We used some integral representations of payoffs of the form <|4.27[) . We refer to [17], [41] and more re- 
cently [20], for some characterizations of classes of functions which admit this kind of representation. In 
order to apply the results of this paper, we need explicit formulae for the complex measure II in some example 
of contingent claims. 



The first example is the European Call option H = (St—K) + . We have two representations of the form (|4.27p 
which result from the following lemma. 

Lemma 4.25. Let K > 0, the European Call option H = (St — K)+ has two representations of the 
form (glFTp.- 

1. For arbitrary R > 1, s > 0, we have 





□ 



4.6.1 Call 




(4.44) 



2. For arbitrary < R < 1, s > 0, we have 




(4.45) 
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4.6.2 Put 

Lemma 4.26. Let K > 0, the European Put option H = (K — St)+ gives for an arbitrary R < 0, s > 

5 The solution to the minimization problem 
5.1 Mean- Variance Hedging 

FS decomposition will help to provide the solution to the global minimization problem. Next theorem deals 
with the case where the underlying process is a PII. 



Theorem 5.1. Let X = (X t )te[o,T] be a process with independent increments with log- characteristic , 
^t- Let H = f{Xr) where f is of the form (3.20)) . We suppose that the PLL, X , satisfies Assumptions^ 0, 
andUH Then, the variance-optimal capital V$ and the variance- optimal hedging strategy ip, solution of the 
minimization problem . are given by 



V = Ho , (5.1) 

and the implicit expression 

X /"* 
<Pt = Zt + -^(Ht- -V - VsdSs) , (5.2) 
nt- Jo 

where the processes (H t ),(£,t) and (At) are defined by 

H t = J v H(u) t ^du) , 6 = j /^dl^Q^ H ^U{du) and X t = , (5.3) 

and 

H ( u)t=e v(u,T)- V (u,t)+* T (u)-* t (u) el uX t - ^ ^„ )t ) =i jf ^*^ d ^(„)_^( )^ . (5.4) 

The optimal initial capital is unique. The optimal hedging strategy pt(oj) is unique up to some (P(du>) ® dt)- 
null set. 

Proof. Since K is deterministic, the optimality follows from Theorem l3.34t Theorem 12 . 251 and Corollarv l2.26l 



Uniqueness follows from Theorem 12.241 □ 

Next theorem deals with the case where the payoff to hedge is given as a bilateral Laplace transform of 
the exponential of a PII. It is an extension of Theorem 3.3 of [31] to PII with no stationary increments. 

Theorem 5.2. Let X = {X t )t^[o.T} be a process with independent increments with cumulant generating 
function k. Let H = f(e XT ) where f is of the form \4-.27\) . We assume the validity of Assumptions® Q 
The variance- optimal capital Vq and the variance- optimal hedging strategy tp, solution of the minimization 
problem |H2]), are given by 

V Q = Ho (5.5) 
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and the implicit expression 

(p s dS s ) , (5.6) 



<pt=tt+4 L ( H *-- v o- 1 ■ 

°t- Jo 



where the processes {H t ), (£t) and (At) are defined by 



7(M) := dpt ^ Z ' 1 ' ) with Pt(z,y) = K t (z + y) - K t (z) - K t (y) , (5.7) 
dp t 

r)(z, dt) := Kdt(z) - 7(2, t)K dt {l) , (5.8) 

A t := — , (5.9) 

dp t 



H t := / etf ^ z > ds ) S*Il{dz) , (5.10) 



6 := / jiz^^^S^Uidz) . (5.11) 



Tfte optimal initial capital is unique. The optimal hedging strategy ipt(w) is unique up to some (P(du>) (g> dt)- 
null set. 

Remark 5.3. The mean variance tradeoff process can be expressed as follows, see fr4.12\ ): 

f* dK u (l) 
K t = / — Kdu(l) ■ 

Jo dp u 

Proof of Theorem ] 5. 2i Since K is deterministic, the optimality follows from Theorem l4.23l Theorem l2.25l 
and Corollary 12.261 Uniqueness follows from Theorem 12.241 

□ 

5.2 The quadratic error 

Again, pdt denotes the measure «<jt(2) — 2«di(l)- Let V, ip and H appearing in Theorem 15.21 The quantity 
E[(Vb + Gt(<p) — H) 2 ] will be called the variance of the hedging error. 

Theorem 5.4. Under the assumptions of Theorem \5.2l the variance of the hedging error equals 

Jo ■= Q J J (y,z)IL(dy)Il{dz) 



where 

p y+z pT 

and 



s o +Z So /%. *. ^)e Ktfe+z)+Q(y,z '* ) ^ : y, z G sup P n 
: otherwise. 



a(y,z,t) := n(z,T) - V (z,t) - (r](y,T) - r,(y,t)) - I [*%Ql)dp B , 



dp s 



P(y,z,t) ■= Pt(y,z)~ l(z,s)p ds (y,l) 

Jo 
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Remark 5.5. We have 

a(y, z, t) = (r](z, T) - r,(z, t)) - (r,(y, T) - r,(y, t)) - (K T - K t ) , 
where K is the MVT process. 

Proof. The quadratic error can be calculated using Corollary 12.291 and Corollary l2.26l It gives 



3 



( e X p{-(K T -K s )}d{L) t 
Jo 



(5.12) 



where L is the remainder martingale in the FS decomposition of H. We proceed now to the evaluation of 
(L). 

Using (|4.22p . (|4.23j) . Remark l2~4l and the bilinearity of the covariation give 

Re((L(y),L(z))) = ±({L(y)M*)) + (W),M*))) 

= \ ((L{y) + L(z), L(y) + L(z)) - (L(y),L{y)) - (l(z),L{z))) , 

and 

(L{y)+L{z),L{y)+L{z)) < (L(y) + L(z), L(y) + L(z)) + - L(z), L(y) - L{z)) , 

= 2(L(y),I(^))+2(L(z),I(I)) . 
(|4.4.3p in the proof of Theorem 14.231 and considerations above allow to prove that 

sup \Re((L(y),L(z)))\ <oo. 

y,z£l+iR 

Similarly we can bound Im((L(y), L(z)) t ), writing 

Im((L(y),L(z)}) = \ ((i(y) - L(z), L(y) - L(z)) - (l( V ),L^)) (l(z),L{z))) , 
so that we obtain 



Im ((L(y),L(z))) < (L(y), L{y)) + {L(z),L(z) 
and 

sup \Im({L(y),L(z)))\ < oo . 

y,z£l+iR 

Therefore 

J J (L(y),L(z)) t U(dy)IL(dz) 

is a well-defined, continuous, predictable, with bounded variation complex-valued process. 
We recall that L t = J L(z) t H(dz) so 

L\ = J J L(y) t L(z) t n(dy)U(dz) . 

An application of Fubini's theorem yields that 

J J (L(y),L(z)) t Tl(dy)Tl(dz) , 
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is a martingale. This implies 

(L,L) t = I I (L(y),L(z)) t Il(dy)IL(dz) 



by definition of oblique bracket. It remains to evaluate (L(y), L(z)) for y,z £ supp(n). 
We know by Proposition 14. 1 II that for all y, z,y + z e D, 

(M(y),M(z)) t = f St z Pdu(y,z) . 



Using the same terminology of Proposition 14.171 H4.21[l says (L(z), M) t = and <|4.20[) imply 

(L(y),L(z)) t = f i e J7«'*Mv«) i j(M(y) ) M(«)>,- /' £{z) e e£ (M, M (y)) , 



e^^ z ' du)+r, ^' du ^d(M(y),M(z)) a - / 
o Jo 

t r t 



e^^ z ' du)+r, ^' du ^d(M(y),M(z)) s - / ~f(z, s)e^ Mz ' du)+ ^ v ' du » S^d (M, M(y)) s 

Jo 



t rt 



Jo 



Hence, 



(L{y),L{z)) t = / ^ z ^ + "^ du »^ z /3( 2 /,z,ds) 
Jo 

We come back to (|5 . 1 2[) . Recalling Remark 15.31 we have 



Jo 



Since EfS'fJ 1 " 2 ] = Sg +z e Kt( ' !/+z - ) , an application of Fubini's theorem yields 

E [ f T e -l*T-K*)d (L(y), L(z))\ = E [ C e<*^ S y t +* (3{y , z, dt) 



o 



(5.13) 



o 



which equals Jo(y,z). 

Another application of Fubini's theorem gives 



T e-^~ K ^d{L,L) t = [ [ ( T e~^ K ^d(L(y),L(z)) t U(dy)n(dz) , 



t 

JCJCJO 

hence 



T 

-(Kr-K t ) d ( L)L ) 



e 



o 



2 



c JC 



e 



-(K T -K, 



) d(L(y),L(z)) t 



o 



U(dy)IL(dz) 



f f J (y,z)IL(dy)IL(dz) . 
Jc Jc 
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Corollaries 12.291 and 12.261 imply that the left-hand side of the previous equation provides the variance of the 
hedging error. □ 



5.3 The exponential Levy case 

In this section, we specify rapidly the results concerning FS decomposition and the minimization problem 
when (X t ) is a Levy process (A t ). Using the fact that (A t ) is a process with independent stationary increments 
it is not difficult to show that 

«t(z) =tn A (z) , (5.14) 

where n A (z) — Ki(z), k a : D — > C. Since for every z G D, 1 1— ► n t (z) has bounded variation then X = A is a 
semimartingale and Proposition 13.161 implies that (t,z) i— > Kt{z) is continuous. 
We make the following hypothesis. 

Assumption 9. 1. 2 G D ; 

2. k a (2)-2k a (1) ^0. 

Remark 5.6. 1. p d t = (k a (2) - 2k a (1)) di ; 

A = ^yj^iy ^W /- an 2/ 1 G [0, T], 2; G £) ; so D =T>, 

3. Assumptions® and\^are verified. 
Again we denote the process S as 

5"t = s exp(X 4 ) = s exp(A t ) . 

It remains to verify Assumption [8] which of course depends on the contingent claim. 

Example 5.7. 1. H = (St — K)+- We choose the second representation for the call. So, for < R < 1, 

I = supp{n)nR= {R,l}- 

In this case Assumption^ 1 becomes I = [R, R + 1] C D. This is always satisfied since D D [0, 2] and 
it is convex. Assumption\^2 is always verified because I is compact and k a is continuous. 

2. H = (K - S T )+- We recall that R < and so 

Io = supp(U) ni = {R}. 

In this case, Assumption^ 1 , gives again I = [2R, I] C D. Since [0, 2] is always included in D, we need 
to suppose here that 2R (which is a negative value) belongs to D. 

This is not a restriction provided that D contains some negative values since we have the degree of 
freedom for choosing R. 

In this subsection, we reobtain results obtain in [31]. From Proposition 14.131 we obtain the following. 
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Corollary 5.8. Under Assumption^ the process (St) can be written as 

S t = M t + A t , 

where 

A t = K(l) [ S u -du and (M, M) t = (k(2) — 2k(1)) f S 2 u _du . 
Jo Jo 

The mean-variance tradeoff process equals 

ft K (i)2 
K t = j r l d{ M,M) u = K{2) ^ Ml) t. (5.15) 

From Theorem 14.231 and Theorem I5.2[ we obtain the following result. 

Theorem 5.9. We suppose the validity of Assumption^ We consider an option H of the type y.ffTp , The 
following properties hold true. 

1. The FS decomposition is given by Ht = Hq + J a T CtdSt + Lt where 

• H t = J H(z) t U(dz) with H{z) t = exp(7 ? A (z)(T - t))Sf and z e I , t e [0,T] ; 

• & = / Z(z) t IL(dz) with £{z) t = 7 A (z) exp(r] A (z){T - t^Sfl 1 andzel,te [0,T] ; 

• L t = Ht — Hq — J Q £ u dS u . 

Moreover, for z 6 suppll, 

A k(z + 1) - n(z) - k(1) 
7 [Z> k(2)-2k(1) 

• ?] A (z) = k(z)-k(1) 7 a (z). 

According to the notations of Lemma \4-17\ we have 

rj(z,t)=ri A (z)t, 1 (z,t) = 1 A (z). 

2. The solution of the minimization problem is given by a pair (Vo, V?) where 

V = H Q and + [Ht- -Vo- G t - ((f)) with A = 



S t - y yr " k(2)-2k(1) 

Remark 5.10. Lemma \2. 14\ implies that is the linear space of predictable processes v such that~E ( J Q T vfSf_dt 
oo. 

Remark 5.11. We come back to the examples introduced in Remark \3.21\ In all the three cases, Assumption 
[H is verified if 2 e D. This is happens in the following situations: 

1. always in the Poisson case; 

2. if A = X is a NIG process and if 2 < a — /3\ 

3. if A = X is a VG process and if2<-{3 + V/3 2 + 2a . 
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Remark 5.12. If X is a Poisson process with parameter A > then the quadratic error is zero. In fact, the 
quantities 

k a (z) = A(exp(z) - 1)) 

Pt(y, z) = At(cxp(y) - l)(cxp(z) - 1) 
cxp(z) — 1 



7(*,*) 



e- 1 



imply that (3(y, z, t) = for every y, z G C, t G [0, T]. 

Therefore Ja(y, z, t) = 0. In particular all the options of type {4-27\) are perfectly hedgeable. 

5.4 Exponential of a Wiener integral driven by a Levy process 

Let A be a Levy process. The cumulant function of A t equals K A (z) = tK A (z) for n A = k a : D\ — > C. We 
formulate the following hypothesis: 

Assumption 10. 1. There is r > such that r G Da- 

2. A has no deterministic increments. 

Remark 5.13. According to Lemma \\~3\ for every 7 > 0, such that 7 € D, 

k a (2 7 ) - 2k a ( 7 ) > . (5.16) 

We consider the PII process X t = f*l s dA s where I : [0,T] -> [e,r/2] is a (deterministic continuous) 
function and e, r > such that 2e < r. 

Remark 5.14. 1. Lemma \3.2J\ says that D contains D E r := {x G R| ex, ^ G -Da} +«R, and n t {z) = 
J Q K A (zl s )ds . 

2. Pt = fv(K A (2l s )-2 K A (l s ))ds; 

3. 2 G D ; X is a PII semimartingale since t > Kt(2) has bounded variation, see Lemma \3.1J\ 
I 1 6 D e<r since 0, r G D\. 

Proposition 5.15. Assumptions^ and\^ are verified. Moreover D e r C V. 

Proof. 1. Using Lemma l4~3l Assumption [6] is verified if we show that 1 1— > pt(l) = Kt(2) — 2«t(l) is strictly 
increasing. Now 

Kt(2) - 2«t(l) = / (k a (2^ s ) — 2k a (Z s )) (is . 
Jo 

Inequality l|5.16p and Lemma [431 imply that Vs G [0, T] 

k a (21 s ) - 2k a (/ s ) > . 
In fact, A has no deterministic increments. This shows Assumption [6l 
2. For z G D £)T , by Remark \5. 141 1. we have 

sup xe[e)r] |« A (a;z)| 



dn t {z) 




K A {zl t ) 


Apt 




K A (2l t ) - 2k a {U) 



< 



infa; g [ e >r /2] (k a (2x) — 2k a {x)) 
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Previuous supremum and infimum exist since x i— ► k a (zx) is continuous and it attains a maximum and 
a minimum on a compact interval. So, D e ^ r C V and Assumption [7] is verified because of Remark T5. 151 
4. 

□ 

Remark 5.16. 1. Point 1. of Assumption^ is also verified if we show that I C D e>r ; in fact D e<r C V 
and I Q U (I + 1) C I. 

2. From previous proof it follows that 

dn t {z) n A (zl t ) 



dp t K A (2l t )-2K A {l t ) ■ 

3. Since I is compact and 1 1— > dK £p* * s continuous, point 2. of Assumption^ would be verified again for 
all cases provided that I C D £ r . 

It remains to verify Assumption [8] for the same class of options as in previous subsections. The only point 
to establish will be to show 

TX 

Ic{x\sx,— eD A }. (5.17) 

Example 5.17. 1. H = (St — Similarly to the case where X is a Levy process, we take the second 

representation of the European Call. In this case I = [R,R+1] and 1|5.17|1 is verified. 

2. H = (K - S T )+- Again, here R < 0, I = [2R, R+l}. 

We only have to require that D A contains some negative values, which is the case for the three examples 
introduced at Section \3^4\ Selecting R in a proper way, (|5. 17ft is fulfilled. 

We provide now the solution to the minimization problem under Assumption [TUl . By Theorem 15. 2\ we 
have 

K A (ls) 



X(s) = 

j{z,s) -- 



k a (21 s ) - 2k a {Is) ' 

K A ({Z+1)1 S )-K A (ZI S )-K A (1 S ) 

n A (2l s ) - 2k a (I s ) 



V (z, s) = n A (zl s ) - kA{21 ^ kA{1s) (k A ((z + 1)1.) - K A (zl s ) n A (l s )) 

hence 

V {z, s) = K A izl s ) - A(s) (k A ((z + 1)1.) - n A izl s ) - K A (Z S )) , 
We obtain finally the optimal hedging 

At /"* 

ft = & + o — (Ht- -Vo - (fsdSs) , 
Jt- Jo 

where the processes (£?(),(£() are defined by 

H t = I e^ v{z ' ds) S^U(dz) , 
Jc 

&= [ 7 0M)e^' )(z < ds) s t z r 1 n(dz) . 
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5.5 A toy example 

Let (Wt) be a standard Brownian motion, we consider X t = W^m, where ip : M.+ — > K+ is a strictly 
increasing function, including the pathological case where ip t = a.e. We have 

E[e zXt ] = E[e zW ^>] = e Kt(z) = , 

so that 

z 2 

= y^(*) , «t(2) - 2«t(l) = ip(t) and «t(z + 1) - K t {z) - «t(l) = 2^0) • 

So 

(M, M) t = j S 2 s _^{ds) and A t = J ^-d (M, M) s = J \s s -^{ds) , 
and the MVT process verifies 

K t = J^_d { M,M) s = p-Ws)= 1 - m . 

All the conditions to apply Theorem 15.21 are satisfied so the function j(z, t) is equal to the Radon-Nykodim 
derivative of Kt(z + 1) — Kt(z) — K t (l) with respect to n t (2) — 2/t t (l), so 

j{z,t) = z, r,{z,t) = ^-(z 2 ~z) and A(t) = A = 1 . 

Hence we can compute the variance-optimal hedging strategy (/? and the variance-optimal initial capital Vo 
in this case 



i r* 

4&t- Jo 



and 



£t = j^{z,t)e£^M S ^n{dz) =f zexp |^y^(*(T) - *(t))| 5^^^) 

Remark 5.18. Calculating (3(y, z, t) of the quadratic error section, we find /3 = 0. Therefore here also the 
quadratic error is zero. This confirms the fact that the market is complete, at least for the considered class 
of options. 



6 Application to Electricity 

6.1 Hedging electricity derivatives with forward contacts 

Electricity markets are composed by the Spot market setting prices for each delivery hour of the next day 
and the forward or futures market setting prices for more distant delivery periods. For simplicity, we will 
assume that interest rates are deterministic and zero so that futures prices are equivalent to forward prices. 
Forward prices given by the market correspond to a fixed price of one MWh of electricity for delivery in a 
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given future period, typically a month, a quarter or a year. Hence, the corresponding term contracts are in 
fact swaps (i.e. forward contracts with delivery over a period) but are improperly named forward. However, 
the strong assumption that there are tradable forward contracts for all future time points Td > is usual 
and will be assumed here. 

Because of non-storability of electricity, no dynamic hedging strategy can be performed on the spot market. 
Hedging instruments for electricity derivatives are then futures or forward contracts. The value of a forward 
contract offering the fixed price F^ d at time for delivery of lMWh at time Td is by definition of the 
forward price, £Q' Td = 0. Indeed, there is no cost to enter at time the forward contract with the current 
market forward price Fq a . Then, the value of the same forward contract S°' Td at time t G [0,T<j] is deduced 
by an argument of Absence of (static) Arbitrage as S^' Td — e~ r ( T ~ t \F^' d — F^ d ). Hence, the dynamic of 
the hedging instrument (St' Td )o<t<T d is directly related (for deterministic interest rates) to the dynamic of 
forward prices {F^ d )o<t<T d - Consequently, in the sequel we will focus on the dynamic of forward prices. 

6.2 Electricity price models for pricing and hedging application 

Observing market data, one can notice two main stylised features of electricity spot and forward prices: 

• Volatility term structure of forward prices: the volatility increases when the time to maturity decreases; 

• Non-Gaussianity of log-returns: log-returns can be considered as Gaussian for long-term contracts but 
they are clearly leptokurtic for short-term contratcs with huge spikes on the Spot market. 

Hence, a challenge is to be able to describe with a single model, both the spikes on the short term and the 
volatility term structure of the forward curve. One reasonable attempt to do so is to consider the exponential 
Levy factor model, proposed by Benth and Benth [9], or [15]. The forward price given at time t for delivery 
at time Td > t, denoted Fj d is then modeled by a p-factors model, such that 



• (m t d )o< t <T d is a real deterministic trend; 

• For any k = 1,- ■ -p, (X^ Td ) < t < Td is such that X^ Td = f Q <T k e- Xk{Td ~ s) dA k s , where A = (A 1 , ■ • • ,A P ) 
is a Levy process on M. d , with E[A^'] = and FarfA^] = 1; 

• (Tfc > , Afe > , are called respectively the volatilities and the mean-reverting rates. 

Hence, forward prices are given as exponentials of PII with non- stationary increments. Then, the spot 
model is derived by setting Sr d = F^ d and reduces to the exponential of a sum of possibly non-Gaussian 
Ornstein-Uhlenbeck processes. In practice, we consider the case of a one or a two factors model (p = 1 or 
2), where the first factor X 1 is a non-Gaussian PII and the second factor X 2 is a Brownian motion with 
&i o"2 . Notice that this kind of model was originally developed and studied in details for interest rates 
in [41], as an extension of the Heath- Jarrow-Morton model where the Brownian motion has been replaced 
by a general Levy process. Recent contributions in the subject are [22|[43]. 

Of course, this modeling procedure (|6. 18|1 . implies incompleteness of the market. Hence, if we aim at pricing 
and hedging a European call on a forward with maturity T < Td, it won't be possible, in general, to hedge 
perfectly the payoff (F^ d — K) + with a hedging portfolio of forward contracts. Then, a natural approach 
could consist in looking for the variance optimal price and hedging portfolio. In this framework, the results 



p 




for all t £ [0, T d ] , where 



(6.18) 
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of Section |4] generalizing the results of Hubalek & al in [31] to the case of non stationary PII can be useful. 
Similarly, some arithmetic models proposed in [6] for electricity prices, consists of replacing the right-hand 
side of l|6.18p by its logarithm. Hence, with this kind of models the results of Section [3761 can also be useful. 



6.3 The non Gaussian two factors model 

To simplify let us forget the upperscript T d denoting the delivery period (since we will consider a fixed 
delivery period). We suppose that the forward price F follows the two factors model 

F t = F exp(m t + X\ + X 2 ) , for all t € [0, T d ] , where (6.19) 

• m is a real deterministic trend starting at 0. It is supposed to be absolutely continuous with respect 
to Lebesgue; 

• Xl = £o- s e- x ^-^dA u , where A is a Levy process on R with A following a Normal Inverse Gaussian 
(NIG) distribution or a Variance Gamma (VG) distribution. Moreover, we will assume that E[Ai] = 
and Vor [Ai] = 1; 

• X 2 = a{W where W is a standard Brownian motion on K; 

• A and W are independent. 

• a s and ai standing respectively for the short-term volatilty and long-term volatility. 



6.4 Verification of the assumptions 



The result below helps to extend Theorem 15.21 to the case where X is a finite sum of independent PII 
semimartingales, each one verifying Assumptions [6l [7] and [8] for a given payoff H — /(soe Xr ). 

Lemma 6.1. Lei AT 1 , A" 2 be two independent PII semimartingales with cumulant generating functions k 1 and 
related domains D l ,T> % ,i = 1,2 characterized in Remark \3.8\ and <|4. 13|) . Let f : C — » C of the form H4.27\ ). 
For X = X 1 + X 2 with related domains D,V and cumulant generating function k, we have the following. 

1. D = D 1 C\D 2 . 

2. v 1 n v 2 c v. 

3. If A" 1 , AT 2 verify Assumptions^ [3 andQ then X has the same property. 

Proof. Since X 1 , X 2 are independent and taking into account Remark l3~8l we obtain 1. and 

Kt(z) = n\{z) + k 2 (z), Vz e D. 



We denote by p\i = 1,2, the reference variance measures defined in Remark 14.71 Clearly p 



dp 1 



< 1. 



dp 1 dp with I, 
If z e V 1 n T> 2 , we can write 



.7(1 



dn t {z) 



dp t 



dp t 



< 



dn 2 t {z)dp 2 



dpf dp t 



< 




dpt 



d Pt, 
—dp 

dp t 



p 2 and 
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This concludes the proof of V 1 n V 2 C V and therefore of the of Point 2. 

Finally Point 3. follows then by inspection. □ 

With the two factors model, the forward price F is then given as the exponential of a PII, X, such that 
for all t G [0,T d ], 

X t = mt + Xl + X 2 = mt + <J S I e- x{Td - u) dk u + a t W t ■ (6.20) 

Jo 

For this model, we formulate the following assumption. 
Assumption 11. 1. 2a s £ D\. 

2. If oi = 0, we require A not to have deterministic increments. 

3. We define e = a s e- XTd , r = 2a s . 

f : C ^ C is of the type ggg j fulfilling (|5~T7|) . 

Proposition 6.2. 1. The cumulant generating function of X defined by 116.20)) . K : [0, T<j] x D — > C is 
such that for all z G D ej - := {x € R | xa s G D\} + iM., then for all t G [0, TV], 

« t (z) = 2m, + + f K K {zo s e-^ T *- u y)du . (6.21) 



2 

7n particular for fixed z G -D e , r , * (— * K t(z) is absolutely continuous with respect to Lebesgue measure. 
2. Assumptions^ and[3| are verified. 
Proof. We set X 2 = m + X 2 . We observe that 

£> 2 =X> 2 = C, k 2 (z) = exp(zm t + z 2 af^). 
We recall that A and W are independent so that X 2 and A" 1 are independent. 

A" 1 is a process of the type studied at Section l5~4l According to Proposition 15.15} Remark 1 5 . 1 61 and (|5.17p 
it follows that Assumptions El [7] and [5] are verified for X 1 . 

Both statements 1. and 2. are now a consequence of Lemma KTTl □ 

Remark 6.3. For examples of f fulfilling l|5.17p . we refer to Example \5.17[ 

The solution to the mean- variance problem is provided by Theorem 15.21 

Theorem 6.4. The variance-optimal capital Vq and the variance-optimal hedging strategy ip, solution of the 
minimization problem . are given by 



V = H (6.22) 

and the implicit expression 

A /"* 

<Pt=tt + ^{Ht- -V Q - I v s dS s ) : (6.23) 



St- Jo 

where the processes (Ht),(£t) and (At) are defined as follows: 

,-\(T d -t) 

za 2 + k a ((z + 1)2) - n A (zz) - k a (z) 



z t : = a s e 
7CM): = 



ri(z,t) 



a 2 + k a (2z) - 2k a (z) 

z 2 a 2 a 2 
zm t + —i- + k a (zJ) - 7 (z, t)(m t + y + n A (z)) 



dt 
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erf + k a {2T)- 2k a (T)' 
6 = / 7(^*)e^" ( *' ,ia) 5fr 1 n((fa) 



TTie optimal initial capital is unique. The optimal hedging strategy pt(ui) is unique up to some (P(du>) ® dt)- 
null set. 

Remark 6.5. Previous formulae are practically exploitable numerically. The last condition to be checked is 

2a s £ D A . (6.24) 
In our classical examples, this is always verified. 

1. Ai is a Normal Inverse Gaussian random variable. If a s < then \6.2$ is verified. 



2. Ai is a Variance Gamma random variable then \6.2^\ is verified, if for instance a s < l3+ ^ J ^fL . 

7 Simulations 

7.1 Exponential Levy 

We consider the problem of pricing a European call, with payoff (St — K) + , where the underlying process 
S is given as the exponential of a NIG Levy process i.e. for all t £ [0, T], 

St = e Xt , where AT is a Levy process with X\ ~ NIG(a, /?, 8, /i) . 

The time unit is the year and the interest rate is zero in all our simulations. The initial value of the 
underlying is So = 100 Euros. The maturity of the option is T — 0.25 i.e. three months from now. Five 
different sets of parameters for the NIG distribution have been considered, going from the case of almost 
Gaussian returns corresponding to standard equities, to the case of highly non Gaussian returns. The 
standard set of parameters is estimated on the Month-ahead base forward prices of the French Power market 
in 2007: 

a = 38.46 , [3 = -3.85 , 5 = 6.40 , n = 0.64 . (7.25) 

Those parameters imply a zero mean, a standard deviation of 41%, a skewness (measuring the asymmetry) 
of —0.02 and an excess kurtosis (measuring the fatness of the tails) of 0.01. The other sets of parameters 
are obtained by multiplying parameter a by a coefficient C, (j3, (5, fi) being such that the first three moments 
are unchanged. Note that when C grows to infinity the tails of the NIG distribution get closer to the tails of 
the Gaussian distribution. For instance, Table Q] shows how the excess kurtosis (which is zero for a Gaussian 
distribution) is modified with the five values of C chosen in our simulations. 

We have compared on simulations the Variance Optimal strategy (VO) using the real NIG incomplete 
market model with the real values of parameters to the Black-Scholes strategy (BS) assuming Gaussian 
returns with the real values of mean and variance. Of course, the VO strategy is by definition theoritically 
optimal in continuous time, w.r.t. the quadratic norm. However, both strategies are implemented in discrete 
time, hence the performances observed in our simulations are spoiled w.r.t. the theoritical continuous 
rebalancing framework. 
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Coefficient 


C = 0.08 


C = 0.14 


C = 0.2 


C = l 


C = 2 


a 


3.08 


5.38 


7.69 


38.46 


76.92 


Excess kurtosis 


1.87 


0.61 


0.30 


0.01 


4. 10~ 3 



Figure 1: Excess kurtosis of X\ for different values of a, (13,8, y) insuring the same three first moments. 
7.1.1 Strike impact on the pricing value and the hedging ratio 

Figure [2] shows the Initial Capital (on the left graph) and the initial hedge ratio (on the right graph) 
produced by the VO and the BS strategies as functions of the strike, for three different sets of parameters 
C = 0.08 , C — 1 , C — 2. We consider N — 12 trading dates, which corresponds to operational practices 
on electricity markets, for an option expirying in three months. One can observe that BS results are very 
similar to VO results for C > 1 which corresponds to almost Gaussian returns. However, for small values 
of C, for C = 0.08, corresponding to highly non Gaussian returns, BS approach under-estimates out-of-the- 
money options and over-estimates at-the-money options. For instance, on Figure [3l one can observe that 
for K = 99 Euros the Black-Scholes Initial Capital (ICbs) represents 122% of the variance optimal Initial 
Capital (ICvo), while for K = 150 it represents only 57% of the variance optimal price. Moreover, the 
hedging strategy differs sensibly for C = 0.08, while it is quite similar to BS's ratio for C > 1. 




Strike price Strike price 



Figure 2: Initial Capital (on the left) and hedge ratio (on the right) w.r.t. the strike, for C = 0.08 , C =1 , (7 = 2. 
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K = 99 


K = 150 


ICyo 


50.08 


7.11 


0.40 


IC BS (vs ICvo) 


50.00 (99.56%) 


8.65 (121.73%) 


0.23 (57.30%) 



Figure 3: Initial Capital of VO pricing (ICvo) vs Initial Capital of BS pricing (ICbs) for C = 0.08. 
7.1.2 Hedging error and number of trading dates 

Figure [4] considers the hedging error (the difference between the terminal value of the hedging portfolio and 
the payoff) as a function of the number of trading dates, for a strike K = 99 Euros (at the money) and for 
five different sets of parameters C described on Figure [TJ The bias (on the left graph) and standard deviation 
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(on the right graph) of the hedging error have been estimated by Monte Carlo method on 5000 runs. Note 
that we could have used the formula stated in Theorem 15.41 to compute the variance of the error, but this 
would have give us the limiting error which does not take into account the additional error due to the finite 
number of trading dates. 

In terms of standard deviation, the VO strategy seems to outperform sensibly the BS strategy, for small 
values of C. For instance, one can observe on Figured! for C = 0.08 that the VO strategy allows to reduce 
10% of the standard deviation of the error. As expected, one can observe that the VO error converges to the 
BS error when C increases. This is due to the convergence of NIG log-returns to Gaussian log-returns when 
C increases (recall that the simulated log-returns are almost symmetric) . One can distinguish two sources of 
incompleteness, the rebalancing error due to the dicrete rebalancing strategy and the intrinsic error due to 
the model incompleteness. On Figure |4j the hedging error (both for BS and VO) decreases with the number 
of trading dates and seems to converge to a limiting error corresponding to the intrinsic error. For C = 1 
and for a small number of trading dates N < 5, the rebalancing error represents the most part of the hedging 
error, then it seems to vanish over N = 30 trading dates, where the intrinsic error is predominant. For small 
values of C < 0.2, even for small numbers of trading dates, the intrinsic error seems to be predominant. For 
C < 0.2 and N > 12 trading dates, it seems useless to increase the number of trading dates. Moreover, 
one can observe that for a small number of trading dates N < 12 and for large values of C > 1, BS seems 
to outperform the VO strategy, in terms of standard deviation. This can be interpreted as a consequence 
of the central limit theorem. Indeed, when the time between two trading dates increases the corresponding 
increments of the Levy process converge to a Gaussian variable. Hence, the model error comitted by the BS 
approach decreases when the number of trading dates decreases. 

In term of bias, the over-estimation of at-the-money options (observed for C = 0.08, on Figures [2j [3]) 
seems to induce a positive bias for the BS error (see Figure|4]), whereas the Bias of the VO error is negligeable 
(as expected from the theory). However, one can observe on Figure that the difference between VO and 
BS bias error is smaller than the difference between the Initial Capitals, therefore one can conclude that, in 
our simulations, the BS hedging strategy induces more losses in average than the VO strategy. 

However, to be more relevant in our analysis, we have compared on FigureO the performances of the BS 
hedging portfolio with the VO hedging portfolio starting with the same Initial Capital as the BS hedging 
portfolio. One can observe on Figure [5] that this approach allows to reduce the standard deviation of the 
VO hedging error (increasing the bias and of course the global quadratic error w.r.t. the VO strategy with 
optimal Initial capital). 

It is interesting to notice that, in terms of skewness and kurtosis, the VO strategy seems to outperform 
sensibly the BS strategy for small values of C. Figure [6] shows that for C = 0.08, the skewness of the BS 
hedging error is strongly negative (3 times greater than the VO error using the same Initial Capital) and 
the kurtosis is high (14 times greater than the VO error). Hence, in our simulations, BS strategy seems to 
imply more extreme losses than the VO strategy. 

In conclusion, the VO approach provides initial capital and hedging strategies which are not significantly 
different from the BS approach except for log-returns with high excess kurtosis (with small values of pa- 
rameter a in the NIG case). Similarly, we can observe (though the figures are not reported here) the same 
behaviour w.r.t. to the asymmetry of the distribution: the VO approach allows to outperform significantly 
the BS approach for strongly asymmetric log-returns (with high (absolute) values of parameter (3 in the NIG 
case). On the other hand, in more standard cases, the VO strategy seems to be comparable with the BS 
strategy in terms of quadratic error and to have the significant and unexpected advantage to limit extreme 
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Figure 4: Hedging Error w.r.t. the number of trading dates for different values of C and for K = 99 Euros (Bias, on 
the left and standard deviation, on the right). 
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Figure 5: Variance optimal hedging error vs Black-Scholes hedging error for different values of C and for K = 99 
Euros (averaged values for different numbers of trading dates). 
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Figure 6: Empirical moments of the hedging error for C = 0.08, TV = 12 and K — 99 Euros (averaged values for 
different number of trading dates). 



7.2 Exponential PII 

We consider the problem of hedging and pricing a European call on an electricity forward, with a maturity 
T = 0.25 of three month. The maturity is equal to the delivery date of the forward contract T = T c i. 
As stated in Section [BJ the natural hedging instrument is the corresponding forward contract with value 
s o = e - r (T-t)( F T _ F Tj for all t g p )T ] ) w here F T = F is supposed to follow the NIG one factor model: 

F t = e Xt , where X t = [ a s e~ x{T ~ u) dA u where A is a NIG process with Ai - NIG(a, 0, S, /z) . 

Jo 
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Figure 7: Hedging Error of BS strategy v.s. the VO strategy with the same initial capital as BS w.r.t. the number of 
trading dates for different values of C and for K = 99 Euros (Bias, on the left and standard deviation, on the right). 

The standard set of parameters (C = 1) for the distribution of Ai is estimated on the same data as in the 
previous section (Month-ahead base forward prices of the French Power market in 2007): 

a = 15.81 , (3 = -1.581 , S = 15.57 , \i = 1.56 . 

Those parameters correspond to a standard and centered NIG distribution with a skewness of —0.019. The 
estimated annual short-term volatility and mean-reverting rate are a s = 57.47% and A = 3. The other sets 
of parameters considered in simulations are obtained by multiplying parameter a by a coefficient C, 5, \x 
being such that the first three moments are unchanged). Table Q] shows how the excess kurtosis is modified 
with C. 



Coefficient 


C = 0.08 


C = \ 
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1.26 
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Excess kurtosis 


1.87 


0.013 



Figure 8: Excess kurtosis of Ai for different values of a (/3, <5, fi) insuring the same three first moments 

Figure [9] shows the Bias and Standard deviation of the hedging error as a function of the number of 
trading dates estimated by Monte Calo method on 5000 runs. The results are comparable to those obtained 
in the case of the Levy process, on Figure [9J However, one can notice that the BS strategy does no more 
outperform the VO strategy for small numbers of trading dates as observed in the Levy case. This is due to 
the fact that X t is no more a sum of i.i.d. variables. 
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Figure 9: Hedging Error w.r.t. the number of trading dates for C = 0.08 and (7=1, for K = 99 Euros (Bias, on 
the left and standard deviation, on the right). 
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Figure 10: Empirical moments of the hedging error for C = 0.08, TV = 10 and K = 99 Euros. 
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